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Elastrerty, Preevelectricity aud Crystal Lattice Dy nance s 
R.U. Mindlix 


Lepart ment of Civil Engineering Columba University, Wee York 


Abstrat — A review of some receuf developments in the area between the 
dynamical theory of cry stel lattices, th tthe harmonce AYyroxima kro aud the 
4 


classecal linear theories of ela streta and preeedectricity . 


(. Lntroduction 


The differential equations of the Classical theory ot the auyiso- 


tropic, elastic continuum are the loug wave, low frequence htt’ of 


the tinite Ar tf ference C4 uattous of the Aynamecal theory af crystal 
lattices of mass farticles . Conse ucatly, Solatrous of the difference 
equations converge to the corre sponding solutions of rhe dh Heveu ta ( 
equations as pertods aud wave lengths (WCVEGSE . Conversely, the two Q- 
verge as periods aud wave leagrh s drurcrish or, 14 the ease ot egulib- 
rium, as dimeuscons diminish. The discrepancy ihvepen greater be tusen 
classical preeselectricity atd mederu theentes of the cry shal lattice of 
polaritable atoms because these two co not eveu converme at the long 
wave, fou Frequency lieth. Tn both cases, rhe detects are due wo insuf- 
ficient con sideration out the couftimuyu theories, et the 5 fratery o* crys. 
als aud the mtevacticus between atoms or molecules - 

During the yust dozen years, there hus been 4 revival et 
activita amed at ex fending the range ef the classicas theurtes ef 
Aastirty aud wero electreuky ti account for more aspects ot Struc. 
ture ans interatomic teractions . Developments have jareg ressed aoug 
numerous paths | Ww cluding : C1) rediscoveries of the arly terms aud 
the while of Cauchy’s Ci] wafiacte series représenfatrou chau elastic 


solid with a periodic structure — eguv alent te Fer king ito acount 


all of the qrachents of strain cn the potentiat encryg density 5 (2) form- 


TR TE MSE TATE. GRO a eke A 


oS pe retene 


ior BF ey WZ 2 RIP SFO PUT FTE ERATE UTA EUPIA TIEN 
aye : . = PRR el ngs SS a THES a PT TAEDY SEY ESTPUORTENAe BE pF ANA “a, 
at FRE GATE FELT EIES T CE OLE RSE PAIN IR RIEDEL IID ER OLS, 8 f i 
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ulation oF Continuum theories of defermation aud poluri zat on et ery: tals ut #4 
Com peund lattice stmcture ; (3) revival aud extensions OF the Cosseral ej 


Contimuum represtitation of molecular crastals ; C#) augmen tatcon of the 


classecal (Vowt (3]) theory eft pieecelectrecity te include the contribution 


EE of the gradientet the electronic pelarteatien, 

fe 

Mot all of the new developments contribute Ao reduction of Tae 
guy betucéa Conftruum aud lettre theories T they @re ytaged on ¥4€ 

: 4 rounds that the extended range 0 f the continuum theorg mast conform to te 
lattice theory in that range aud thal rac auguedred ami Yagt at: 
: COmmedutée bhsrved er olscrvable Phasical phencmena net ucccunwes fer wy 
rt.- Clarstut thé or y. However, the Following Favorable coucluscows Can 
7 reached : 


Ud) The Cauchy- type He on tes are ace Cpta ble tt ak feast Fae first 


Aud scrond gradvents ot shan are included . Then the addetronal terms in 
q the dis persion relation for plane wuves Can be matched ty the carly terms ‘ 
% ot te lattice Ais persion relation wethouf sacrificing positive oe Kivick- ness 

q of the potentcal Cusugy density; aud the phenomenon of surface cuergy of 

a 


detormation mekes ifs appeariuce for contosymmerric as wellas mon - 


Centresy mm etree crastels (tte latter peg aicins only te stra ano its 


q first gradient}. 

4 (2) The contimuuns eguatcons develused vocently fr Me teuece 

3 crystals ave the covrect, loug wuve (units ot the difference eguations of 

the Mu C/- tae lattice of both muss parteles and polarizublé utcurs 

i and vhey give tHe correck be huvior of the ophtal branches at long 
: Wave leugths , Aatcommuda te sur tac€ enevay of de formumcu (ana polav- 

q (zation) aud, with rhe iiclasion of the magneree Kitld, qld a de> - . 
4 persion relateon @ xh rocking the ong wave orteons of the caryled Ccouste, 

: optical Lue electromaguetec branches . 

4 (3) The eg we toons of the Cosserut continuum are the long 

. wave limet of the ditterene equations otmotiou of monomolecular 
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crystuls and contribute the long wave behavior of the soft optical mode. 
The recent extension of the Cosserat thearg fo include a micro -structure 
that is deformable, as well as rotate ble, may prove to describe the huh- 
er modes. 

(4) The addition of the Cust gradiont ot rhe elechoure polaria- 
tion to the energy density tn classical piexoekctreity su pyle s the 
missing terms reguired to match the long uuve linet of the di tfer- 
ence eguatrons of the lattice of polarreable atoms ant serves ty 
accommedate several observed phenomena otherwise nut iacluded : sur- 
face energ 4 of deformation and polarization , anomalous capacitance ot 
thin dielec thee films, acoustical activity aud, wrth He wiclasron of 
rhe Wnagn ere Field, optical activity . 

in the Following pages, typical examples of these deve loy - 


ments aud results ave deseribed within the Framework of linear con- 
frnuum +heones (linear di Ffevential eg uations } aud lalfice phewtes 
in the harmon a yproxin atoon (linear difference Ca eatecns }, The 
pont of view «sy Frow the side of continuum theories and they are 

ft -ated 14 more general terms than ere the lattice theortes, For t4e 
latter, refereuce mag be made to” Theory of hal tice Dynamecs 14 the Harmon 


Ayproxin eke by Maradudin Moutrofl, Wess and L patova ( Ace demre Fres, 
1471). 


Tn RA ALE DEM ON A NPESER IE LTE FET ot Magee te Cag et 2 en dS ATEN TEA ge SURES, NE ABEL ETN Nes Sa BERET OER RT hanes te ER TALE log Lene ELA TERN De OTSA ES 0 Bd BAS BANE DAI ER ALE Ft OPES hie 
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2. Simple - Cubtc Lattice of Mass Ponts and Classical Elasti:ty 


The fact that the equations of class cecal elas teity ure the long 
wave, (ow Meneses lint of the difference & uations of a lattice of mess 7 
particles, indicates that classical Elastic 1s tnnted te wave lengths aad 
bodily dimeuscons large in Comparison with the distume Yetuen uecarest 


newhber atoms (the “Jattne parameter”), For the actual magactudss of 


Sra MERI Cee NII Be orl vend ates Oe Be nae namin Bs 


DAA Se LEE TET Pah BN RST Pe DEY GENE icant, 


the errors it us necessary to compare detaled solutions cf unaloquus 
d 


pe 


problems in the twe topes oT theory. For thes Purpose s/s conven 


reat to employ, for the lattice So canes those ftrmultarce by Cu 215 


ELS OD OPI 


: 
x 
& 
3 
a 
$ 
we 
} 


Herman ana Welles [4] fer a Simple cubic lattice with nearest ana nepl 


néaregt neighbor central force interactions ano ther novel angular 


SEP PESOS LOS 


Interactoons ue tween three, Successive non: cullincar efems . Simple- cubic 
; is the simplest of aif lattice structuves and the Sars -Hermaa-whius Es - 
9 uations ave the Simplest valid, cubic -lat tice cg uatrons that do aot 
ers a relation among the three elastic constants Ef cube a aaa try 
Pritts ug b ne natural crystals with simple - cuore structure ure kuowa to 

5 exist, the cg uatrons are suituble for « study of the masoret fects 


ot wave length and Slze. 


i, Dif rence F quattons . 


ta rectangular coordinates KE 5st : (,2,>, the atoms ot rhe 


Ste 


‘ 


: ead 
pees sta 


, 


simple.cubic lattice arc Taken fz be ut the points 4 tu, Xoima, &y7 


Biel once 


a Na, where Lun are positive or Acgutive Integers and 4% Is rhe 
b, laftie purameter. The ceutral ferce constavts (€, fore por uml 

a relative ds placement) betveen neavest neighbor utoms and botween nexh 
mearest neighbor atoms are Aesignuted by & ana /, respectively , while ’ 
y 1s the ai haus for< constant between three SteCess IVE MOT vilintar 
atoms. The csles of the three force cousfants arc Depicted schematicull 


7 
by dushed lines tn Fig. L, 


ear 


: 


ETRE eR OE I BT TEE De AY OTIS CEE ER PR Ae FABER UB TE ee cee wane rece 
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With a, t=1,2,3, the rectangular Components of displacement ‘ d 
of the atom at £, wn, Gaeis, Herman and Wallis [4] find three equ “ston 


ot motion Cf the type 


M u hen S+l, mn wubb 2 hen) 


=< Cu, i : 


RelLmetn i-\m-tn kel min R-l,m-l,n 
+ elue Pea NO ea run! 
satan * Rel net fetpm,a-i helm ned ry 
rU ™, ru, 7 tu, mM, +u, am, ~ Sui") 
hiimein t-t mein Rel mein 8-L mela 


«(bey iie' +U, ~ U2 Uy 


Ag lm, nel 


R-\ mm, a-t Sel m,n-i helm. ve! 
+ Uy +g ane -~Uy alles “Uy wees ) 


4 =f a -t & : 
ty (ur. um ares nal tute Au, ee (2.1) 


Where Mis the mass of the atom. the remaining two Cg uutions vt notion uré 
Obed by Syctcal permutation cf subscriyts and Sepe stains. Born doe ven | 
Karrie, tr there + cers ical paper [5], woulé have found he Same tgucethens 

Fite had fate ser 64, 6,6 fo be the present a, ty, 3,0, Brg, respcctiely, 


At free boundaries £-+L, the concittious to bw sumsfied are ; 


tog (ao Mea Sees 
plus be Any u, "(eel mtn * ute, ae 2 eel) m,n . gutters) ; 
ra FS ( Big tee _ gee Z Ops ae hi ~ us be “, Het) | 
£ ar ata urbe : yo * urtemnet : authms ) ‘ 
ry ge tales 2 Us terl meta -ustemeln , y themeta) 
a(usees m, nel uy orm nel ust, nel, Genres) iO: ‘ (2.2) 


1. Reducthen to Equations of Classical Elasticity. ‘ 


To show how (2.1) reduces to the corresponding te ua ton ia 


Classical Clastety | IT 5 convewrsut fo putit na different form lel, Define 


difference cyeraters as follews. 


rs 


gz w - 4 , 4-4 r “ 
TA = Career -u; », Ya, 4; ura (ue a u, ™, "Va, . 
7 t mm, a) 2 ; 

a an 2 a % ae e Cat mn seus " uf ome) Ig ; (2.3) . 
7 

Rim, + + - = 4u,n ; ores 
4,4,4,°"" -b(alah +a,a; +070, +4, oh) as” i 


’ 
cy 
SES ret ees wre te ed Anes hy tae stad hea 


RIP TES FOE CIE EEC I REE BINS PAO MP e Lae OSE eA NS 


Peay 


ie 


cd 
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and analogous defiatons, e.g. ayasens (ag cae uo" \/a, for forward , 


: 
| 


vackwurd , eeu cen wal anal cross differences in the remein ing coovdinate 


' divections Thea, noting that 


Feloreln fe-lLmein £41 metin Rel, nse in Pre) 
u: Fads ha Re a = a*latatat +242 424% U2, 
1 e 
hel mein £-Um-t an k+l mtn £-lomet &@. 4) 
GCS BE enn eal ~ as &a* A Alu. o%" , 


} We Can write (2.1) in the form 


sca fet 


. =legh + Cqg( Ap +4) * $a “GA *(a2 +A Mus mt riciatCy) (4,4, by eA A345 wat (2.3) 


! 
where 


‘ - ro - . 
PHM/a, acy sr Atha, O6,.268, Acqa? <(b42/), Le; 
: j ' 


fan 


Mow, expand the difference cperuters, eo: © ul") i Taylor series ef par tel 
; 
derwwutives 0: , = O/dxK;, , Of continaous disvlacement funchons u:(4,) iTS: 


i 2 mn _ 1 2 2 2 ‘ 
Ai 4; =(i+ a 9, +e ) dru, 
Ban if, 4p ee tye a 2p 
a. A; jak =(j+te o+ga : te) Bay (2.7) 
: ina Ltd gb *) 4.8, 22 
! 
[T¥€ cnly second derivatives arc vetauined, (29) reduces to 
oo 20 2 2 ; 
pu, = Ce,a% + Oyq(22t2s)Ju, + (Cat Guall 2, d24,+ d,d3u5), (2,8) 
which ¢s the 2g uatron ot motion of classccal elastreity for muftrials with rhe 
Couns tants Cus Stes Coe of cube Sy mucter 
1 
iii, Dispersion Relations, 
; The designation “long wave fow eres weneg limit os for (2.8), sfems from 
a compartion ot allspersion relations ( freguency Vd. wave number) For plane t 


waves from (2.5) and (28), Consider for Examuyle longs tadeau | Waves in the 


' X,-Airection . ‘For (2.57) we set 


4, mA ima Lyman. i(gha-wt 
. 2 -0 A= Ao i y 0<ge% rT. (2.9) 


ee fa D8 DS i a es oe Si eel ean 


a LOREEN Red , Lie OCS SES Ser aon > CANIN READ ae eg TSAI SF TAOS 3s WERE BOY IN yc 845 AR PRS, PRES THAR LEER Teen see bee ts oat eet o- 
Soiet § 6 Sees Rg Tee nes < NRE 

ee 

ta 


a 
SAR 
et 
RE as 


rary 
ps 
~ 


¢ 
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Bt “4 


Wh 
es 


Then (2.5) Is satisfied if 


Sits 


ph 

ex a 
‘ 
- 


n rie (2.10) 


Q 
=~ 


For (2.8), we set 


UzsU,70, ur Ae Groot) (2.01) 


and find 


Syebpeesrustee ed bs cea etree tiie: 
AN a REE ee tease es 


w= 6 (2 (2) 


EES EON 
eo 


The two dispersion relations , (2.10) and (2.12). ave illustrated im Fig « 2. dt 


’ 
7 


may be seen tat the lattice dispersion relation (2.10) approaches the continuum 


2 relation (4.12) as the wave length approaches intinite , ce. as Phe wave numyer 

eS 

Es upproackes zero. Simultaneously, the Frequency , (ne monatomic Brava lattic: 

Phat 

FY: ra 2 

ES also ayproaches z2¢ro, 

‘e For more detuléed examinations cf the difécrences $ctucen solutions of 

a lattic€ ana continuum éEquukcons, th is NEC € 5ar4 to consider problems for 

‘eal u 

a 

a bodies uth at least onc Fincte Apmension with whch to CE meyer c rhe luff 

Me 

ES parameter a. 

ie iy, Thickuess-Shear Vibratwns of a Plate (8), 

ay Fee Cag a a Oe ee me 

e In a plate bounded by Leet, consider the Qs Jlaccments 

Bs 

oR ma 4 ” . . 

*3 Pe aed: win = (A cossha +A, singtale’*, O< ga<ir. (2.13) 

%) 

A With these displacements, the First two eguatcons of the type (2.9) are satis fied 

a 

a identically aud the third is satisfied if 

‘4 

ee 

‘a paw = tc,, sin'tsa. (2.14 

4 Upen substituting the displaements (2.13) in the boundary conaitions (é.2), we fad 

; 7 that the first two conditions arc satisfied tdentreulla aug the Hard 1s sabes 

14 fied if 

me 4 
pe gas pr/(2i+l), p= 0,1,2,..¢2bel, Q.is) 4 
RS ; 
$ 


ty a i a % ee eA ee EE Ds TL Sie Fe OP OES EEO EET ST UBER SOE FE ISS UE 
Fes * q ra “ BEE £27 & or yaaa x SS FRR ee ORS BITE ET es Sete ro < a 
* raha 1 ty et Pye os LEE SAE? riser reed SW) eae oe FRR LATS 


soured 
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where @ven and odd p apply to symmetric and antisymmetric modes, respectwely. 


From (2.14) and (415), rhe Frequencias are 


2_( Sel" 5, Pm 
ws 2 (Ss) sim 5 e=9,42,...¢ abel (2.16) 
e 2(2L+1) 


and the displacements | from (2.73) and (2.18), arc 


Laan 


uy’ = (A, cos al + A, sin ere \ it (2.17) 
abel att 


To compare the lattice selutton with FkE Continuum Che, We Consider us 


continuum plate of thickness Zh, where 


Zh = (2t+l)a, (2.18) 


be. Each layer otf atoms 1s replaced by a thickness a of the continuum | 


The frequencies of the thichness-shear mides of sucha plate arc gwen by 


r (12 
e w: SES) ; rr 6,1,%.. (2.19) 
e and the displacements arc 
a U,= (A, cos PO hy gia 2 je , (2.20) 
2h ah 


In the long wave, [Ow et ayproyimation, 


TR terse Foe 


fa-2%, pele, (2.21) 
: in Whteh case (2.16) and (2.17) approach (2.19) ana (2.20), respectiveiy . 
P. Lf the Sréguenci ct the lowest, antisqmme toe , thickness -shuar mode : 
: in the continuum plate, 6.6 pel in (2,19), (s used as @ reference Frequency 4 
fe s 
' “2 
: ww, = = (4) z Tea See\ "2 P (2.22) 
a 2ah\ oe (Zrtija e) 
¢ 
3 the normalized Frequencies for the lattice are 


5, 


SF aie pts” Cevid eesn pew eel Loch 


es 


PTR Ee TE Sei ree Sa ESTE IN nA PRCA Ae Nai Meme RA oo AS OANA a pen Ee Bele RLY OAL, MINS Oe SONI ET RE OOS ASE SME LA ga fd Be ero aoa i... 


% 
a 


Seer eer 


a(2bel) pr 


in—_—— =01 .. &2L4] 22 
wT 2(2b +l} ? r s ne; 4 ¢ 3) 


ud 
Qhem- 
w, 


and these are to be compared with the normalised Frequencies » PAO 2... 
for the continuum . 
The discrepancy between the Frequencies p and 12 and the 
Corresponding departure of tne lattice mode shape frou the sinusoda( form 


Charucteristte of the COUTINUUM, as the order of the mode increases are 


Uustated in Fg. 3 for the case of fifteen layers. It may be seen tat, fer 
the first few modes, where the lattice wade shayes are nearly sin usordaul the 
normal ceed srequescies pand fl ave nearly the same. Insecction of he case 
ped shows that the half wave length Can be as small as five lattice peru 
mefers with ce good representution of the mode shuye by « sinusoid und 
an error in fuse tenes of less thun 2%. However, 45 the mude shupe 
departs from the sinusotdal for ta creasiny orders, the Frequcmees pand 
nt. separaté — by almost F¥oOV%e for the highest Mode | tn thease of 
fit teen layers. As the number of layers Increases | the percent discrep: 
ancy between the Cregweneics of the highest lattice mude and the 
Corres ponding continuum mode approaches SOM -100= 57 Yo. 

It should ve noted +hat, whereas the number sf moses in the 
Coutimuum plate ts unlimcted, the number of modes in the lattice 
plate ts equal to the number of layers. This (s because a wave described 
by N pots con have no more than N-I nodes. 

V. Face- Shear ancl Thicleness- Twist waves tn P tate [8]. 

Waves ina plate, with gisplacement and wave normal af right 
angles fo cach other and parallel to the fuces of the plate, are called 
thickness twist waves except fer the wave of zero order (in wuich rhe 
disslaceneat does uot varg through the thickness Of the plete, which Is 
called « facc- shear wave. In the continuum, these waves arc special 


cuses of Love's {9] transverse waves Ma super Ficcal layer 


SOF I OAR Ae Ve TEE SEDER ET TRON NON HAD DEQ LENEAMEP DT RENT OR EM EN DEEPA TO EG RNG REINS NMI TEER at OTN ANT om Oe POT Sa 


FRE FEO FE ae et. FP! 
; 
3 
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In the lattice plate, the wares ave represented wy displacemcuts of 


the form 


» Us" 2(h cos gha rA,singlale ime? *) : (2.24) 


where OK tac O€nael. This time, the eguatrons of motion cf the type 
5 / 4 rf vid 


(Z.9) require 
; pat w* = 4c,,( sin’ dea + sin’ bya) (2.25) 


aud the boundary Conditions (2.2) agai require (2.15). Thus, the normalised 


: frequenc 1€s are 


w _ 2(dhel 2 pir 2 tg ) 
2 2 = S52] sin eee sin Lte -012...<2L+\, (2.26, 
Xd we ( a(zu+i) ¥ ’ a4 ) 3 r a0 ? o 


PRET Ey Bes 


Sip Bincnetet be 


and the displucements ave ‘ 


wee = (A, cos ae +A, sia ey ee : (2.2() ‘ 
atel 2uri 


Tn the long wave, low Frequency limit (La —k,, Ma X,, «ael) 


DPtnied EEE eee ELIT 


4 tHese become the kuoum results for the special case of Love waves [9]: 


} neal pteleghAl™ , ps 942. (2.28) 
4 u, = CA, cos(prrs,/Z2h) + A, sin (prs, /2nyie 1% a (2.29) 


The real branches of the dispersion yelatcon (2.26) ave illastreted im Fea. 4 
E for a plate fifteen layers thick (L=7). At infinite wave length (4-2) along 
a the plate, the mode reduce te the tuickuess -shear modes illustrated ia Fig .3- 
¥ These variations of desplacement across the theckness cf the plate are muin- 
: tamed for all wavelengths , 2r/y, along the plete. The major differenes . 
between the d's persion relations for the laffice and ie cadens 


a Cach of the finite number 2Ltl, of branches of the lattice despersion 


. relution (2.26) hes a hegh eeg sence cul of F rT) add Tien to a low 


Se LTT A 


it ; 
frequency cut-off (2.23), aad all have the same upper cal of € of wave 
Number 4> la, Ce ware leugth dasa fo 2a; wherras all the infinity of 
veal brauches of the continuum dispersion relation (2.28) are © hyperbolse 


curves extending From low Re geeney, cutoffs Q=12,...,0, at: se 0, +o 


infinte Frequencies and wave numbers — asymptotic +: 


7 1 
N= (2L+l)qa/mr, (2.30) 
‘ 1 
which is the straight line starting from the lower left corner and passing 
Through the upper right cornev of Fig. 4 dnd ise tect, the face shear : 


branch of He continuum dispersion relation . 7ha's, the Sis perscou rélaken 


for the continuum IS agood approximateon to that of the lattice chly v9 
the flower le ft region cf Fog. 4 ; ¢.é. for long wave lengths, w COmpurisan 
with a, both along the plate und across its thitkuess , 7wo er three times 
a ts sat frccenl for es al within aie 


i 
Asin the case of thy Cbinine iin, the lathe tapiaites relation 


has branches Yor real Freguency and Imaginary wave Widuaieey eek ee net 
U t 


show ® 1M Fig 4. 


The problem of face- sheer and thickuess- shear waves in a plate 


has also been solved for the foce-ceatered cubic lattice { Brady £10]) aud fer 


the body centered culce lattice (Gong fu). 
vi. Torsional Eguilibrean of ao Re ctanguler Bar, 


The bars bounded by free faces af 4: tL and m=2M ands ia eguilib- 


rium under a twist about the apis of ¥y with angle of tast T per unt length. 


By analogy with the St. Venant solution of The Couations of classical 


elasticity for the analogous problem Cra} it is assumed that 


= Ta* (tm +A sin £6 sinh me) O<O<T, 


Peat oe 
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rik 


bi 
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eee 
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With these al plecementsy “the hae two egustrins of the type (2.7) wrth He 


left hand side zero, ave satistied saeeiaally and the ive! eguation ts sate fied 


i 


wih if : 
cosh yp = 2-cos6. ; : (2.32) 
a ' } 
The conditions of the ae (2. 2), for free boundaries reduce, wn View 
of (2. 3), to bs : i 
| ; t | Hi 
22 (ustrimn _ usher) ylerne! - Oa =0 on = +L, : (2.33) 
Gated, nae, & Ayn n= 
2 2(ue ‘s “Uy a "\ ous mune <a =0 on m=tM., (2.34) 
iu i : . ' 
Dpon substituting (2.31) in (2.33), we find ; 
i) | | ! , 
q 6-8 mere Vw /(aust), pr Lash. (2.35') 
’ i 1 { 
s Heme, the turd of (2.31) may he wriften as a finite series + . 
! . . ' i 
; ke, 2 vay é : ; 
i us" = Ta*( km + ZiAy sin £6, sinhg, ) , (2.36) ‘ 


: Substitution of (2.36) in the boundary condi tions (2.34) yields 


e=b ' 


: ; ! 
z Ap sin 46, sinh 24, cosh(M+ £)¢, Sel L562 jy ks ; (2.37) 


C.@. @ set of L simultaneous,  ltacar, alge bree equates on rhe coef fiir ‘5s Ae 
% 

The aah of eg uations Can be Solved explicitly for the Ap by « method ana/—' 
: Syhs? te. that for Serene Fourier coefficients | Multiply beth sides of 


(2. $1) by anit and Sum over L from kz 1 es -L. Mow tisd, 


’ 


fobs SLO, — bcos (L +40; 
BSINLG, = ——n ze — EE é.3 en 
: i ae 4sin $6, 2 sing 5. ere! : 
1 ‘| Also, employing (2:35), we find ; a j ’ ‘ 
4B ' , i ‘ * 
' a me 0,4 ¢ 


E, sin £6, su 46, = 


' 
j : : 


(2.39) 


' 
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hd 2b sin $6, cos (L+3)6, - sinLt, as 
(2n+; ) sinh’ $4, cosh( M+ 3) 4, 


Substitution of (aro) in (2.36) Completes the soluhion tor the Wer pie, function 
eee, The results ave shown in ig. S i which the ratios of Wategers a7 
the lattue porwts give us’""/ta* at those ports . Fug . S 13 to be Compared 
wit Fig. So -ihic’ depicts typrce! contours of the Waryins fuuchons for Square 
and long recfanyu (ur sectrons from the St. Venant solution. It may be seen thet, 
with one side of the cross section vcstneted to three atoms (L =1), the warp — 
ing bears little yesemblance ty that found by St. Venant » namely there « 
no warping ot the squere section (Ht, m=) aud , ia eack Succeedine, 
section (L=1, M=Z 3, 4), the dis placements at a bng side (45 #1 me 
0 ti, 2, sap do not have a marimum and qminmum between the centr 
aud the ends as they have in the continuum solution - However iT tke 
Smaller dimension 1s increased to fice atoms (L: 2, M=2,3,..) rhe Sissirmilaritg 
disappears For the sguare (L=2, M=2), it may be seen thar the cross scetron 
*s divided into cight sectors, with alternating signs of displacementts 
instead of the usual four sectors for long rectangular sections, This ts 
preasely the result found by St. Venant for +he continuum, as ilushuted 
‘Mn Fig. 6. for the rectangular SeCtons L>2,M=3 and L22,M-#, rhe 
displacements along « long side réackh a maximum ana a uaimem Alar 
the ends, while along a short side, the dis placemcats vary mane ton cally. 
Again, these results are the samc as the corresponding ones in St. Veaunt’s 
solution, as llustrated in Fig. 6, Thus, the cross section need hare onl 
as muna as Five atoms along the shorter side for the displecemcats 


te have the mejor gualitative proyerties found th the continuum salatoon, 
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3. Strain Gradient Theories 


Extensions of classical elasticity To account for crystal strc turl 


began in 1851 weth Cauchy's Cll wnfratte series represen tation of an is oFropie ; 


material wirk a periods structure. Although Cauchy did not carry bieabx 


TAY EIST HS 


very far, if is now Known that hus theory corresponds ty uu intiaste series 


A 


pth 


Cx pansion of diféerence operators to a labre theor terms of ait fereate| 


4 


wperat ors; or, eq uivalentig , augmentation of classical elas fierFy throug 4 the 


Were Bs 


tacluscon of all he gradients of strat in the potental energy 7 de mss ty | hitte 


attention was pard to the work for many yoars— possibly due to the compler: ry 
intrduced by Cuuchy*s consideration of non-centrosymmebre isotrepy } port bla 
due ty anerror of emission in the constitutive 4 teeta for the tymmetrne yor] 
ot tte stress * possibly duc fy todhunter's assessment: “This consists merely of | 


generalities, and (5 «ppartatl, of no importance “ (12, Vol. pi dred. 


i. Rotation Gradient, 
elit acl bnctecsk ED 


Laterest in the gradient tyye of theory was sfimulated, 1 1960, Xv 
Aero and Muvshinskce Cig], Grote [5], Rajagopal [1b] aud Truesdell and Tou UT] 


who took nto account rhe first graaent of the rotation: 


where Sek 1s the unit alternating teusor: 


+t, tyk = (23,231, 312, 


aK -U) tgk = 32, 203, 132, (3.2) 
0, otherwise. 
A connection with crys fal lat fice theory was noted by Krumhans! [/8), ie : 


1963, who showed, by expansion of rhe differcues in the Cxpresseon for the 
potential energy of « general Cravas lattice, that the second erder terms, 
in the vesulfing serves ot clerivatives, contain the rotation gradient 


The novel feafures accompanying the rofation gradient are the. 


nine Wine n= aden rin dbs 


Appeerances of an antisymme fac partot the stress Cindletermmiate J, & couple stress and % 
maternal constant with rhe dimenson of length ; bul, except for q tenuous comecton 
with molecu/ar crystals (see Article 5), no relation to ery shel lethce theory ayoetrs 
#0 have been established beyoud Kremhansl’s brief remark, 
it Furst Strain Grachent, 

The rotation gradient has erght independent components whch are eight 


of the erghteen Components of the first gradient of the strain. The additional 


ten Commponen fs are 


40,9. UR & OF Ou, + 0,0; u;). (3-3) 


Altecnatively ; in terms of the strain ; 


Si, = 5(,u;+2;u:), (3.4) 


all erghteen Components may be expressed as S54 5 out perhaps the sumplest 


form is the second gradient of the displacement : 


Stik = Q: O: fel . (2.97) 


The fundanmenuta/ eguatrons for the elastre material in Which account ts 
teeken of the full first gradient of the strain were given by Toupin [19] in 1962. 


In the linear case, the potentialenergy density may ¥C expressed in the form 


We Cui Sijnt ECijka Sey See 7 E Sip k dn Sif Samant ijk dan Sez Seton 1 (3.6) 


where the Cy, are material constants. The mest interesting feature of (3.6) 15 
the linear term Cig Siig wheek correspends to a self- eq uili brating int fre ( 
stress and gives rise to a@ surface energy of deformation * an observed 
physical phenomenon not contained in classical elasticity . 

The surface energy of deformation 13 @ part of the energy associated. 
with the Separatian along @ surface. Partot the siparation cnergy is ¢he bond 


energy — that required to break the atomic bonds across the surface while the 


AEA YEE LEU ELIE ERAT EDAD DDE ie. 
€ x 


Ree tag oa TIINS 
PETTITT OMS URE POSEY 
; s . UE TA TINGE VES oe 
near RATA TL I EE oe EES SET ae ee ae 


to 


relative pivitions of the atoms in cach of the two resulting portions are held 


fired -- Say, by fietitvous forces, The release of these forces 15 Accompanied Oy % 


deformation localized near the Surteces, and an associated (negatice)cnergg the 


" 


high as 46% of the bound energy [20], The extreme locel(gahtou of de formaticu 


surface energy of de formation — estimates of the maga tude of which range us 


al the surface was first observed in (96/ by Germer Macleae aud Horf mun [2] by 
| means of low energy electron diffraction measurements afnicks surfaces. 1 he 
found that the displacemenT of the sapertiisal (ayer of aloms foward rhe mittrior 
was Five times as large as that of the nest feycr. Toupin aad Garis [22] on 
1963, found mathematical solutions for a similar surface dleformattou within rhe 
framework of Toupin ‘s strain gradien! equations und ulss cbibted rhe currespondeuce 
with a oue-Aimcnsconal, mona toms. lattice with nearestand net nearest neg kbor 
interacttoas, 

The first Strain gradient theory , however, suffers from two deficicacns, 
{ First of af, it con fans a sur tact energy of defermurcua only for non céntrosy meme 
ric materials . Thes 05 evidenced by rhe third rank tensor cocfficient of the linear 
term «4 (3.6); whereas the elasticity of ceatrosyume tric mulerials Can ve char- 
acterized only by tensors of ven rank. Secondly, rhe disperseen relatron for 
plan€ waves Can match tut for a lather beyoud He lincar range only it the 
req uurement ot & positive Strain entrgy density (s Gaandoned - a5 Shown in the 
3 Following scchion. Beth of these detects Cis aypcar if rhe Scion graaénf of the 
: Straw 1s incladed. 
iit. Second Stran Gradient C23]. 

Lf rhe strain and its first two gradients are taken ity acecual the 


potential euceyy density tor a centrosy mmetre malénral hus the form 


q = 2 Cth Si; Ske td Fojkhin Site Skmn 4 Béjk Onn pog Disa Smngg * Pajkkunn i; Shtoun 


q * Bima Pujne » (3.7) 


i 
© 
t 
f 
4 


where, in edition to (3+) and (3.5), 
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IT 
Sijkd 7 0; ;O¢Ug ‘ (3.8) 
This enerqu aud the usual kiaetic energy deus: ty 4 pd, oe lead ty the stress eg catrons 


of motion 


BT 9 — 2625 Tyg + 5 Taj = PMs, (6.9) 


where the stresses Ti; lie aud Tove are given by 


& 


Tyj # PW/05;; = ine See + Fijietnn Skeonn » 
Tk? ow/ OS rig = Kise dm kun 9 (3.10) 


Tae? ow/a Sie - Vojk Lan Sun f Pik Liunpa Sunpe t Pike ? 
and te three, vector boundary tachens | 


q@ a . -~ 
tee Tg BT, 4 RM Teigdd + Lely (Tyg OTeigel] + be Ly (me Tejee) 


~ Lt i, n; nel Dane) ning Tesco | ’ 


(3.01) 
Q . . ae - T . e - 
ty = ncn, (T;:4 Tipe) r nls (nT eike) + Ly (nen; Ucige) 
Q) _— 
cn new, My Vi jke. 
where 
L,= ne (oj -nznw) Opn; - (6. -nen,}2; (3.12) 


Cr willbe chserved that the Sstramn euseay density (3.7) has no term 


linear th Sy. . This 1s because, just as im Classical elasterty , the malerial 


ej 
tonfigurateon, to whith TW ts referred) Can be chosen so thal sucha term 
does not appear, Also tieve is no term in (3.7) linear in Sask as tts mareria/ 
Coetficcent would be a ‘tensor of odd rank which Caunel ces! Cer ceutrasym 
metrre matirials. This leaves the calg linear term to ig the cuc ia Sikes 
which can be shown to produce a sur face enerby of deformation as fullews 
From an ung framed body under no external forces, reniore a portion, 


of volume V, bounded by a surface S. The energy i sucha portion ,6n 


taal ilorian under wo external forces, «5 


wefway, (543) 


DARE, Sy Sncde Bamae nS aey Fh o OPES 
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a 
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or, from (9.7) anq (3.12), 
~ 4 o . 
Wire j (Te; Seq * Tae Sige * Tejas Seite jav es (; Seine Sync EM. (5-4) 


The sCeonal tarigral in (3.1¢) enters beruuse, as may be Ken iin (3.10), af Corteapinds 
to a constant tirm Poses “w Te see whertas the PE wnat ce ty pacts a AL aud Nad 


Vary tA pre portion to the struca ang its gradients ~ the first tn tegra lin (s.14) 


aan bo tauvfermed fu 
Hae Levu. et (nda; + t"(n-d.)- 4; Jd5, (403) 
which vanishes since the surfuce Ss free of fraction. Thies leaves 
WEL Be ye deje AV = 2f5 Me BEine 2 2_ He dS (sv) 


as Qn Cnergy tat yemuns on V al though the body to under no cyternal fortes. 


The Energy per unit arca of the surface 15, From (3.6), 
Wee £5, 22502, Meus | ( 
7 ON Seppe % Me Ae ds: ay) 


To show that this energy 1s localteed Aer the sur fucé, this su flrecenl 
ty consider the cese of the half-sypuce x, 30 with xK,-0 free of truction 
Assume 


U,zU(4,), U,eUy~ 0. (3-18) 
Then the sftss € ua hous of cgurlioriaen , from (5.9), reduce ev 
Dig. Bhan ee tin 4, (315) 
while the Youndury condi trons, From Cli), beame 


7, - 2, 


te £00 Tuy =0, Tut -9, be 7 thas :0 on X,- 0. (3.20, 


ul ; 


Miso, with the essumed form of aisplacements (3.18), the constetutwc 


4 


equations (3.40) treke the form 


Ty 2 COMe yp Duy, Tt A274, , Tuner pa% re oF uy, (3.21) 
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Ee and the sur face energy per unit area (3.17) becomes, after noting that nN, -—l and 
‘ a Biiks* 90> . 
¥ ° 2.2}, 2 
a zc W=-26 ora), ‘ (3-22) 
as Substitution of (3.27) m (5.19) gives the displacement Cpuation of equilibrium 
we ie 
a 
E (¢ - (K-27) 9) +6 be 1974, = 0, (3.23) 
Bb or 
7 i bt) 52 F 
as (1- £722)(1- GO") Oru, = 0, (3.24) 
BG where 
a q , 2 ; 
aA 2c g = k-2y * [(a-2p)?- 4c) pele: (3.25) 
q Regarding boundary condi tions, iT may ¥e scen that the first of (3.40) 
§ : 1S Gutomme teal y satisfied cf the stress equation of eguili ortam (st) 03 hes 
7 fied. he H maining twe boun dara conditisas, from (8.20) und (3.21), arc 
4 ie 4 ., - 
CG -pj up dtu] nae 0, Cy du, +6234), 9 = ~ Bo. (3.26) 
The solution of (3.2 +), Vantshing af in finite , ts 
. us A, enh +A, ehhh ; (3.27) 
The values of A, aud A, may be determined from thc boundary cong, hom > 
(3.20) and they tle surface encrgy of deferiatron 15, from (5.22), 
W*: “2 fo( A, 8" +A, 4°). (3.26) 
‘os The associated stray iniinishes exponentially into the mftrior with oe cay 


rate conteolle d by the magartade of 4, and G. The ayplarance of material 
constants 4, with the dimecuasion of lengts, (s Fy pical at the gradunt a 
theories and yields a length scale characteristic ot the structure ot the 
material. Sucka seale 1s absent inclasacal elashety . An estimate of th 
order of magaitudle of the length seazle may be obtained From He 


- electron adiftractcon measurements at mckel surtuccs by Germer, Mache 
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and Hartman C21). Ther fia ding tat the displacement et the super trial 
layer of atoms ts five times as large as Hat of the next layer would lead, 
on the ussumyp tion et a Simple €x ponentre/ decay @ MC to an f of about 
five etghths ot the dis tence ten adjaccut layers oT atoms. /hus, th 
defor mution 1s qa frcen ¢ only at the first few layers. 

| Upon departure intoy rhe mterior, the displacement nced net muimtara 
thr sr ne 51g 4, The conditions tor positive strain ene rag density do not 


include relutrens between hk «nel 4 or be tuten Auaay, Mewe, trom (3.29) the 


q muy be complex, Corts ponding he an oserl(uturg détay ot displacement taty 


+h tuterier : 


Ty obtain the ‘Lorre sponding solufton for 4 lathe Consider u“ 
Strigl€ row of partcles distrbuted along the X-aKrs, af unct destene 
upact, with teractions as far us third nergh bers anid including sc lf -€54r/- 
[rating mitial forces, Wow, separate the laltice betucen the parteles ef 
az V0 and n-l and consider the scmc-inbiurte lafhre nod. Ihe separ - 
ation ts c f4fectcd bu the addition of forces Be, P, on particles af 
O,',2, respectively, eT ual and oppesite to the resultauts e( rhe iutial 
forces exerted en thuse particles b4 the particles af -[, ~2, -3. Since 


he iutial Gorcés are self-equilrbrating , 


Po+P, +P, =O. (3.29) 


With force constants £4 4,,%3 for the mnteruchens between 
First second und third neva hbors, respectively, and with Uy, the dis 
placement of the nth par ticle) the egucligriam of that parhcle 1S 


CX pre osc d by 
G 


=} 
Z Kilns ~2uU, +uUyi)20, n33 (4.30) 
qs 


ten 
Bd Wag -2uy + U2.:) + dy (Us-U2) = Pe. n=2, 


A,(U, -2U, +4) + d(Uy-4,) + d,(U4-4,)2R, Meal, (3.31) 


d,(u,- uy, ) + dy (4.-U,) #44(U5- Uo) ca A n=0, 


: £8 4 
_ e We adopt the novation s Bug = Une, — Uy 5 é 
a ‘ § 
: E Aun = Una 2U4n thy » Ody = Une -tUne TOUn 4g tUy.2 . (3.32) : 
: F : ; 
: a Then . ati th 4 
-_s Byrd tthe + 9d, B22 -og- by 5 Asedy, G.33) i 

3 : ; 3 

i the genera/ equation of eguilibriam (3.30) becomes os 

¢ . 

iw 

J (6.- 6, 4* +4,A*)au, = 0, (3.34) 

L 

. 

Ye 

i (1-A2a7)\(1- AL) Aru, =0, (335) 

where 

2. 2 /2 . 
22,4; =},% ( 8, ++ 8,85) ? t=h2. (3.36) 


The soluhon of (3.35), vanishing utiafiaita, ‘$ 


-nd 7 
u,=A,e +A,e"% (3.37) 
provided that 
f 6; = arccosh(1+1/24,), c242. (3.38) 
3 y As for the boundary conditions, there are three for the two Consteats 
q A, and At, just asin the case of the continuum solution. Fallewinc Gazcs anc 


Wallis [24], we sum the three equations (3.31) to obtam 
(1A? ai(I- Ag A) Au, = 9, (3. 39) 


. which is sates tied identreally by te solution (3.31) with (3.38), Thus, asin 
the case of the continuum, te general Condition of equilibria d's puses 
of one of the threo boundary conditions. The remaining tuo conditions 
Serve te determine /t, and Ag, 


The Sum lars ta in form between the differential equations , Pounder 
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Conditions and Solutions for the continuum and the lattice cs ty be 


remurked cend the correspondence of the continuum expressions to the 
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22 
long wave wave behaviur inthe case ef the lattice may be verified oy noting 
that 
Aru, = OU, + de Opa, + 3 Of, wee 
Atu, = dtu, i : dfu, boos (3.40) 
Beds Of Fe 
and, in (3.38), 


9—> ifr, , (3.40) 


It is also necessary to compare the long wave behavior of the dicyer- 
ston velations for the continuum and the lattice. The one-dimensional Cqua 
v 


tions ct motion, for the twocases, ures for the continuum , from (3.23), 
Ces (a-2y)2} + oO* {dru = 4, (3.4.2) 
and, for the lattie, from (3.34), 
( P,~ 6247 +6;a*)A‘u, =m4, . G.43) 


Lnserfing the wave Forms 


U,= Ae’ Wn-ot) Oe 


3 n“n , 


(3 4+) 
vespectively, we find, tu the fourth powcr of the wuve numver, 
gw = cs" Fe (K- 2p) §* te yj (3.45) 


Mw* = 8, 6°- (3 ar0)0- +7 . (4.46 ) 
T€ onl y the strain audits first gradient were considered, (245 ,weuld be 
pw: Cereus 


and k ts reg ired to be positive by he Condifion of positive ce femsteness 
ot the encryy density . Heme, the first gradient theory) the grea 
veloorty, Aw/dg, mast increas€ as the wave Number MCrtases trom 


zeru, but, in the lat fice the opposite 1s genera 4 trae, as vila strated 
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in Fig 7, In the second gradient theory, o 1s replaced by &-Zy as the cocf- 
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fiereat of the second termin the dispersion relation * and the conditions of 


LE posttivé definrteness pace no res ttetions on the relation between of ana py, So 
pri? 7) | 


ie tlut the group veloeties rey ether increase or deerease wit Wi cree sting were 
4 num br. 
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“4 
4, Diatomie Crystals and Compound Continua 


L€ the laftice structure of « crystal (5 & Stmple one of the 

Bravass type, with one atom yer cell there is no distinction between tHe 
long wave limit cond the long wave low freq uence limit ot rhe furcke dif- 
ference equations of motion or the dispersion relutton, ttowe ver tf the 
luttice has a “basts L 6.0. & group ef two or more afous pe cell, there 
are Wigh trea weaey “optical branmke> as well as /ow fregueneg acoushe 
brunchéy ut rhe (ong ware lomet as illustrated in Fig .E, In thet case, clas- 
sical elasticify accouuts for the long wave porttons of only the Gcoushe 
pranchis. The uygproprate extension [6] of classical elasticity to accom- 
modate the long wave behavior oF ig fclyareme crystal Is Suggested b, 
the fact rhat @ latthce witha bests can be ve solved into? two or more inter~ 
penetrating Sravars lattices. 

The simplest lathce with a sus 1s the Na -/ tape which has part 


icles af the yrorats (la,ma,na/ of @ simple-cubic latter ; bul the particles 


/ 
ut Lemen Even have mass mM, sayy aud thoscet £tmen odd have mcs MA, 
as represtuted in Fig. 9 by curcles and dots, respectively» J timay be seen 
thut the lattice can be resolred into two, saterwenetrating, facé-Cénterid, 
Cubce laftues: cue with particles of mass M and rhe other with parkicles 
of muss M. The wttrachions 9c fucen particles may be teken ty be Stucler 
te thuse for the Gazis -Herman- Wellis lattice described in /rtcle 2.i. the 
force constant vetween né€arest neighbor (unlike ) particles (Ss designated by 
ah. The Force constants between next ucaresl ncighbor Cuke} particles ar 
assuméd toe Ye diffcrent: 3 and 2 for masses Mand A respectively, Ihe 
angular ferce constant, {. '3 taken to 9c the same whether a parkicle . 
of mass ™ or M is at the aper . Then Garis and Wallis C25] tind , for the 

particles at Lemen even, three equations ot the ty ye (2.1) with Mand @ 


re placed by ™ and (: and, for the partecles at Lemen odd, three gus tons 


ot the Type (2.0) with Mand 3 replaced by rm] and me Whew the Bis place meats 


See eee 


Sune rs e RRS PEELE TOTES Ep aes Sai 
for se 


Oe, 
& 
rh a: ‘ 
RNAI ARNO, OE Nn BOA ae mete aren rmnenrtuNRAnmaNE Latina ee NOREEN we gai cee nena eva 6 


as 


4m, v C4 emOz end -wt) : 
u; . = A, 2 > J ; (*.1) 
ry : 2 \ 
with Vzland 2 for Lemen even and odd, respectively, ave substituted in the 

‘ : : J ; i 
Sly oiffereuce eguations of motion, there results the dispersion relation 


toa 
dy dpi diy dy, O 

los ry ' 

| dy daz day 6 dy, 0 : 
soe 4 


dx dyz dyp 9 0 day 


2 2 za = 0. (4.2): 
dy, 0 ) dy, du dis 
e a 2 : 
0 dz O da da das : 
a 2 a 
C 0 d3y dy dy da, 
: ' 
where 
v vy» 2 v ae 
db, =Me -2(d+By)- +6 (2~ cos 6; & cos), ee . | 
d; = ~4( Bey) sin 6,508; , LS, (4,3) 
é di; = 2ncos6: + By 2 2 cos 6, : ; 


im which & jandk range over 42,3 and Vv over 2. 
I¢ we desig nate u, and @: as the displacements ot particles with 
mass rm) and Mh, respectively, and write , ) 
ecu | 2 UG, +m 8,418, -wt) 
i 
fer all Ama, and adoyt He notation (2.3) for finite atFfevence operators, 


» . ms 7 1 
the sawe ais person re lation as (4:2) results from three eg ue tions e t, the 


type 
is Mayet” = 2(a+ Bae = ee) tZa 70(24; +A,+A 5) a, een 
haat had | Ayman : 
+#(beylot (Aa, a" +A, asus") 
; fymyn 


- pat (a, a, + AY A A) ee, 


. 


Aman 


tata af per at lay +ay)d, 


and three more obtemed by interchange of Superseripts laud 2. Ths 


form Le] is more convenceut for 4 851ng ts continuuar fiuits suce the 
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(4.5) 
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Mt : necessity for distinguishing between Lemen @ven and odd is dispensed with, 


1 At long Ware Leng th aad low Frequency us a. SU; , Sey . Then 
I 


¢ { 


€ mplo ying the Expansions (2.7), discarding derivatives higher than rhe second 

: and adding correspond iny members of the two sets of tarte cg uations of 

ro motion we recover rhe Cquations of Hesieatelesticrt of rh tape (2.8) 
with 5 aaa , , 

a p= (mrm)/2a? (46): 
] \ ' 

| aud stiffness constants’ given by (2.6) with 

4 ; 5 ’ . 
cq gf hy  eetieaie: ey 
| . 

q At in weve lengrh s but not Mecessarily low Pregieneg (4.57) becomes ; 
a swith’ m= 2a? py 

i ab yg 2 ake bpd, 4) 5 (228 + DE eB), 

| | + 2(bey\( 0,2, u, +0, 0,45) 

| +ha 2 a, + ay (a 1%) di, “t ' (4.8) 
4 Thus, for the long ware ayprérimation , rhere ave fhree eg uu Rons ct the 
| t4 pe (4.8) aud three more chtained by iatercnang ¢ of suyersertyts laud 4. 
3 : ‘The sry equations «yield the long. weve limits of the dispersion relation (4.2) 
for bath the acoustic and optica! brauvches. 

A continuum theory that produces cguatcons of rhe same form as 
) ; (4.8) may ve conspructéed ds Follows (6), We consider two inte gene hebing 
J rs ideati fre dg ba superscripts, land Z, representing the tur? fuce = , 
| centered cupric sub-latties th the Mall s,fructure. The potential! Enema 
7 density is taken as ai Guud rutee Functron of the strains cf The two 
4 Continay and therirélative diy placement? qud pretation « [er rhe crastal 
q ; cla $5 mom to which the face epukees bute add Visit lattes Oe Lory, 
: : this is, with #512 and A=b2,' ; > 

E weet(5,.-s - 555 tealTuguy ec ars we ¢ 22,¢ eS Sin (4.9) 

Rite i ee co i a a's : } 


al hd ws Alb earns The KONE le ra ey 0 Te 
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4 A oh, Shs She eee ee, " i 4 2 : eee +1 Ae ne Me Ae RAL 


sevefeannjabseyincii ele benwenat Ren ee 
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where 7 


ee fey ee 4 ¥ % 4 
u =(u-u,), Pa > 5 (2.4, +0-U,), a; = 3 (d;u" -d:u0), 


. KA _ AK _ {KA eA kA KA 4.10) 
“eke = Cejka = (cits ~ 2e4e) Gijeat Miz 5; bee + Cea (Si, 8g * Die dix) 
aud the ie are ung tf alliadues are alike and vera vtherwire. 
With Kinetsc euergy deusits 
ae . ee oR ; 
re 2 = (ieee (lt) 


and the potenhal energy density (4.9), Homilten's vartational principle for inde 
¢ 


‘ ont ra _ 
pénteat variations bu und bu, quields the six Fuler cguutions 


FY 3 -(-1*La?*(a, -u,) . OF 2,0; ~9; a) 
tt, A 


) 
A KA A A 
+ ELce 0,0; U, + Cag O;(2; 4, +d, a.) | : (4-12) 


These become the stx egus tions of the type 8) 1f we set 


(cK KK « KK aR 
Aly alan = 26, Aen 2 Otly, cv =O, 
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(4.13) 
2 2 a+ 
Aly F- ay mit ae. ac, 7B4, a = (a+ By)/a?. 


thas, the theory ofa compeund Continuum, reyriseare by the 
potential and kinetic euergy densitifs (4.9) and (4.11), produces Couations 
ot motion which ure the (ong ware lumel of the Caxts-Wallis €inite 
dof ference equations ofan Mad type lattice of mass purncles , 
eg manu oc noted that the linear term in (4.9) «5 the cory g 
. density of « se1F - CG urli brated inctial stress which pruduces u ser face energy 
ot deformation analogous to that ia the second sham gradient equa tens 
described im Article itl, but with the continuum réprésentathon of only 


nearest aud ACKE nearest neighbor tu teractions , 


“ae 
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5. Moleculay Crystals and Cosserat Continua 


Crystels with identical groups of closely packed atoms, Form to molecules 
situated at the points ot alattie are fermed moleculer ery stuls tf rhe 
Mtramolecular forces art large n Com yurisen with ehe intermolecular forces 
(265. Ln a monomoleculer crystal, there 13 On molecule ateack point ofa 
Bravers lave, The as placements Of the mass cenvers of the moleecles corres- 
pord to the dishlacements of the aloms of u momatomic fg f free : and the 
Assocceted €4 uations of motion reduce ty those of classecu/ ela shes hy a 
the long wart 4 pprokimukcon , Hou ver, the rofuton aug strain of rhe 
molecules themselves may alsy be taken into aceount, La rhe lone, wure 


approximation, these mo roms , whirk may be \termed the micro-rofulton and 


a He mocro-stram of a mitve-strugture (molecule ) | Cun be represented 64 ay 


4 Asymmerree tensor of the Stcond rank) say 3; ‘ The AA sym met ree url 

3 of the tensor : = 
eS Yeajr = 2 (46; =e) (3-1) 

4 describes the wmicru-votation aud the fymmetric yert: 

¥ 

§ a 

ig : -2t . ‘ > 

EE tj) F 2 CF. re) 5.2) 


describes the strain os both of whiel ia d be differeal trom tke rotation 


ond strain calculated Erou the dts placements Of the muss centers, 
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Ihe nitro-rotution was faken inty account i a Continuum theory 
E oy Ev and F. Cosseruf [27] in 1909. Their eguu fens were revived and 
applhed to twe- dimensional problems by Schae ter [28] 1 1962. Su bseguent 
Iq) rte Cosserat theory was termed “micropolar clasticity 7 by Eringen [29] . 
Comparisons ef various forms of the theory and gentral Sduttons of the Cf wetrons 
by Eringen [29] , Aero and Kavshinsk ee £30, 31], Meuber 132] and Mindtin [335 were 
given by Comm [34,35] along with a more concise gencral sedation. As ts shown 


below, the Cosserut theora ts the long Ware appr ixidicatean toa lattice theory 


iS Soo Bi el est as ee VE gee a ag asl tag 
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of monomolecular crystals . A direct correspondence betwean the Cosscrat tcory ad 
the molecular crystal KIV03 hes recently been exhibited by Askar [ 36]. 

The long wave, low Feeey lime oF both the molecular latte and Cosserar 
pheoriés 1s the rotuteon gradient theory meatesed in Article 31. Thes linet rg for ut 
1s the samc as thal obfained from the Cosscrar Theory by con Sr ia sig te 
Micro -votation 0 be the same as the ro fatto Caleulered fron the displacements, 
in fact, tougsa [375 refers to the rokizon gradient theory as the “Cossernt theory 
with constrained rotations: 

Extenscens of the Cosseral theory to tuke wato account the strug 
ot the méro-structwe were made in 1964 by Brigen [385 (imicromaggshee” rheerg ) 
Green und Rivlin [39] (multipolar wheory ) and Mindlin [40] (“micrs structure ” theory). 
tn ther Simplest form, these themes are Charachenred. by the potential 


cnerg4 deusety C4o] 
W-3 Coke 3 SegF oz b ike Lj (ee 2 Fei bean Koike an 
t Diikem ly qeem & Fk ue tT | Senn TG ike Kj Src , (3) 


where Si; ts rhe ordinary stra (macro- straw)». 


Seb (2.u, ae) (64) 


ia rs the relative defarmation (the difference between the displacement gradient 


dud the mtro-de formation): 
fej 7% - Vij SS) 
and Werk is the gradient ot the micro-de formation’ 
iit = Oc Bik > fo.e) 


Lu the centosymmetac, Usotryne case [to], Fi jhbm 88 Frpty , 17 (5:3) ave zcro 


and 


cone = AG de (Sip Sia t es din)» 


seach oe AS Ware Be a tee nh eee ts sla a en ceeela EN tat tweet 2 


eres) 


Cee taht AA EERE Bh Re Cron + 


q 
: 30 
‘ Gejee = 1 Se Sea G2 Seg Ge + he Ue), 
bijke = I, See t br Sin Gig rhs Cig dg » 
: Qijtdmn = &, (i; Ske uci ay 2 Sn Sen) ta, (v,, Sine Ont + Ue Vin Sim) (5 ?) 
tay; Stn Shy t Ae a Ji San t Qs Sok Sinn Sut +b Wy) Sun) 
4 45 GS Sac 1410 Sit Gim Sen 74 (U0 Cem Sin + See Yim Tin) 
; ta, Org On One tag Os¢ Cun Cin rary Ons J: tis, 7. 
4 The associated kincké encray deasity ts 
e e 2 e e 
T= 24, 4; +¢ed ve iy 3 (8) 


where ? 13 the mass density aud d 's « characteristic length dimension ct the 
molecule (half the edge length in the cose of cube). 

jhe Complete setof dispersion rclahens for Plan waves in such 
mettrtal with a defornable miivostucture 3 given in [eo] aqua ckhigts twelve 


branches: the usual three acoustrc branches And jin addi tron, nine cyteg/ 


brauckes, Me. Branches with NO sero Fregucney al Infinite ware length. At- 
low freguencees und long wave lengtss, the ditterenfhal cguurons of motion 
aud he d tspursion relations reduce to those of the Firs! Stren Graces eal 


theev 4 ' 


Iu molectley erystals, the most wteres ting ot the oyftca| branches 


ne 2 


ts the relativelg low Tong asec “Tibrational” or “soft optical” brunch (lade tea 


TRO, fer transverse ro tuhtona( optical in [+o$) taki coupled Wy a transverse 


acoustic branch, the eguatens of metion for this pur of couglee mudes are 
ebtuned by setting tua) = 9 (there by elitr ena fing the shear mitra strain) in tthe 


first and third of Ess. (6-5) in Ceos uth the result 


-~ a = ee 
ky Ue ~Byg % Yoay = PUR » 
= - a - y ie 
Way Pat ky Foray 2 Kis ty = z VF Faq > 


where - - 

ky = t2g, +b, y Kis b,-b 

uw ~ Gt 22 "43 2 73? (5.10) 
Ky3 = 724g Gy Mgt Ley -~ 204g Hig Hy» 


Re 


CS ane I eee FS OE cee ae eT ee 


ee, 


ia wet we a TP caw ers ne ev aS 
‘ 4 
nev nraes a Gv settee char te < 


3! 
Since Yyy= 0 in (8.9), theg ace equation > of @ Cosserat coutimuunt. Setting 


uae me Hee | deg * Bette) su) 


anol eliminating A aud B we ob fac the dispersion rélafion 


ks -~pw* kia 


7 SrA ee Re ee (5:12) 
Ki38 hy + 2k, > God w 


The two branches of (5.12) have ordinates gud slopes af $=0 as follows: 
a ae. ale 
, 20, se] -(2kucke\” 
acoustic branch: ow), a 3 hee ( 2 
(5.03) 
tical b h: skis s des =0) 
soft optical wranch, a pa ee 


it may be vert fied that (5.9) are the loug wave limits of the fate 


di f ference equations of motion of a linear tattice of dumbell molecules 


L2.. Artiule 10). We consider a line of molecules of length 2d, Spacing a, 
muss m and momcutoft raertia I with transverse displacement u, and rota- 


tion Ww, af lattice pont M, as tllasprated ty Fig. 10. The potential energy ot 
the lattice may be taken as 


¥ id a z [34 ,(4,.4 -Un) * bau, ~Uy-dYn) 


- é 434° ( Gnot -$,) taky (Une ~ Un nes “U, -dy,)] ($. (¢) 


where Lissolg are force constants, Com paring with (63) we see tat (6.4) 


(5 the one-dimensronel, fiucte difference analogue for he coutosymmefee 


Seika 74 Feit P42 
Su, = (Deeg; + jut; (Unga, ij 2 2G Fig > Umer “Un S44) A, 6) 


Kijke a Hi" (+ net - tna. 


With fwite Aftereuce operators A’, A” ond 2 again defined as im (4.3), 


the eg uations of mofron ure given by 


rig i ae Ge a ol i lg Ok a i Rl a a 


— og = 2 ‘ 7 
Un (lta, 42d) au, ~ 9 (d,rdg) A Y, = muna 
e oy _ 3 + 2,2 d> oe 32 es 
3 y, ae (agra) Un rad OY, ~ az 42%, > Ip.q 
lasevt ; 
ae Bet lna-wt) yp ieee Bat ena -we) (ey 
In (Silo) and fiad the despersion relation 
2(k, td, r2K4)(Cossu-1) — mut Dla, +dal(cos ka-1 - sit Su ) 
rf ae =0 (6/8) 
hytdg)(cos bu -1 singe) Zd,d*(cosfa-I)~ dd *e Ie 


which (> Mlastrated in Figs ll. The oramates und sloyes of the fee branches 
ar g=€ ure: 


acoustie brunch: wy 47 0, 4 = ete te tera) 


4¢ ¢:0 md, a° 
| a 
: soft optical branch! ol zd a ’ =| =¢0. 
20 L de 
$20 


ff, tn the Fimcte ditherence Cgc trons ot mofton ($.l6) we set 
4 


ak, 2K rk,teg y ahi 2 (kth )dla =a,d [ha , aky, > k3d*, 

E 3 2 tg (5. 20/ 
Re - & 

Ape eee 


aud relaca only the first Term wa te eK pau sion otcuch of the ditFercnce 
operators in series of Afterenfal operators, we Find that (9.06) wedece 
to the di fF feren tral equations of mokou ($9) of the Gssercl continue. 


demlacla , iT we insert (5.20) in the exyresseons (5.9) fer the long ware 


: JS epee 
RoR de ST 


huacts of the two branches of fhe lal tice disyerscon relateor , we fies Me 
Corves ponds ne EX yresstens (5:03) for the Cosserul coutinuune s 


Aw e@iternative formulation of the lattice ae and anu appl - 


cation to KNOs are given ky Askar ee 
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6, Classical Preroelectricity i 


The classical, linear theory of elashe qielectnes (prexoelectricity) is 


expressed in terms of twenty-five vanables [3 41] : 


U; Mechanical dis place ment 
Si; Strain 

Taj Stress 

E; Electr freld 

De — Electruc displacement 

P. Polarization 

pe  Electne potential 


These variables are linked bq twenty -five equations : 


y (a) stress “4 uations of motion 
1 
: Tae +f; = OU; ; (ost) 


ae 


(6) equations of the elechostotic field 


; De 20, (2) 
E,=-%, 53 (6.3) 


sy 4g) oa 
(d) constitutive relations 
€ 
: Ty 3 Cisne See ~ Cua E., (6-5) 
= $ 
Di Sipe dug + Ej, EK (G-4 | 
P= Di-e Es, (GY 


where 


ate a » AEN OA NR 


OE OP CT ae PE Pee CP LT ee ee 


RPI ENA PE ORE FTI LIOR PTI LOE eatye eR ES PAO SLES ATS S&H E 


EAE A aT TEE RON ee AMY, CLE RPV TE EOE AEA ey WamL OE owes then tree sea EBay 8H te 


y 
R 
te 
Doe setae ate TARE 
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Cite = elastic stiffness (E; = constant) : 
Crk = piezrelectric stress constant : 
ey, > permittivity (5,; = constant) 
é, = permittivity of a vacuum 


Be ip pk ecm Larsen ace Pan RE cancel 
~ x 


The equations ot motion muy be réduced to feur by eltusuau toon of 


all tHe varce bles except U, and yi 


Piatt I Sc allt ca 


€ oo 
Coie Abed * Ske Puc * F7 = PU» (6-8) 
$ 
Seej Maju Seg Bay =P (C9) 


A linear version ot a varcational principle due to Toupin [+2] may be 

employed to produce equivalent equatcons oft eguilioriam jas follows, : 
First, separate the eneryy deusity, W, of the dielectne tuto the ‘ 

Stored energy of deformation wad Polarization, WS and aremunder whieh 


| 1s the < nergy Geusity of the Maxwell electne self: Field: 
W-Ww'(5;),P:) tse 4i 4: (6410) 
Then oletine an electnce enthalpy [4] 
H=W-E&,D.. (out) 
Upon substituting (6.t0), (6.7) und (asin (641) we finel 
Hew" (5, , P-F eobies te iP; - (6.12) 


Consider boda occupying a volume V bounded by « surface S 
sepavatiny V from an outer vacuum V'. For such a Sao tem, the verseon 
of Toupin ‘5 principle that leads ty equations eg uivalent to those pre 
ceding , ter the case of equilibria 3 


Fh Aa 3% 


- Sf Hav +f (bu, rEPSP avr ft; bad5 =0, (6-13) 
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for independent Variations oF ue, Pr and. Ta (6.13), V*= VV’ and tts the 


surface traction across S aud E! 1s the external elechec Field. Wow, 


t b H 
wy 65,54 PU OP, ~ 654s Ses + Gz SP, HP See G. 14) : 
4 ‘ : 


Define a stress T;; and an effective local electne fore EF by 


c t N 
T; = aw J Eis _ ow ‘ (6.15) x 
: oS; oP 


Sep £ (Uj; uz 5), (6.16) 


PPR Ae HE Ld ede ae OF Mee oe 


Ty Oe; * 18; Su; 5 > (1,.6u; 


oA ie vj,t 4) » 


€,8e = (¢; 84); - Gis Sy 


P. oY; = (Pp; S¢},; = Pei fy , 5 
by the chan vule of differentiation. Then ; 
BHT; bu; ~(EZ-@;) dP, ~ Cages rF i) 5¢ tl T, dur ew; +?:)4¢)); « (6.17) 


Votwig that Vi=VrV' ano thet u; and ?; donot erst nV", we find, 


after apply sia the divergence theorem, 


“6 Hay = [Uj r(EZ-@ oP + (-€, eet P. bpd 


= [gngl Ta Suse Cesgrr O—l ds, (0.0) 
~ bf HaV=-f, bot SPAV-[ eames seds. (G16) 


Hence, (6.13) becom s 


{ 


ge rs) du; +(e: -@; + £5) 6P, +(-64,, rh, bp] dV -/,6€j hed 


t LG, “nT ey) Su; +h (e.he,1-P,) se | ds=0, (6.19) 


where Te} ¢s the jump in £; across 5, Then the Euler cquations 
3 y) 


(6.20) 


(6-21) 
a and the natural boundary Conds trans 
d by J on S, (6.24) 
4 n(- ey J +P.)20, 


follow from (6.19}. In (6,22), n; (-e,¢ +P) ts the suv face charee. 


the energy clensity ot deformation ana polanzutice is ited i €e 


re. gp AS PD, er -?P 


2 “tyke “ij See is Fee; 95; VK (@+23) 


$0 thet, from (6.15), 


t $ 
-E: ts a Ph + f7 S 
J wk’ jkk PRL > ar 


T. t 
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Equuteons (6.46), (6.20), (@:20) and (6.24) with boundurg con- 
Aitwons (6.24), Constitute a Imear version ct the (quutrons of egusli beens 


of elastee clielectnis in the form given ay Icuyin Le). 


; $ P 
The relations between the now constants aj, Pek aud Cepee ate 
the TKE Standard [43] ences ave found us fllews , 


From (6.3) anol the setond of (6.20) (omit ting Po aid the 


E usual fermalatcen ) ES 2 Es, tteuce, the first of (6.24) becomes F 
4 
4 EAP Fee Shee (6.25) 
y Vow the ratio of IP to e,F ts the ditlectrce suscepti bi ity, Hence ae, 1s 
4 ‘ pope S$ 
4 proportienal to the réccyroce | susceptibility and constant Stra (x3) “ 


$ a -| $ ¢ 2 
Ay ~ €o Kij - (6-26) 
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esata d aie gate 


min ET atER. tg Ste eke EARP, Ay ge 
See EMRE REE ARS RD A EMPRESS AR re Bee BON AIRES ENE ROE TS PLES 


SRSLECTROT AONE aera nee stn paneerneners ance 


ce ne aA, 


| 31 
1. 
a wit (6.26) ; (6.25) becomes 
ae ° 
el 7 
| ’ Es = & Xik Pet Fine Dee - (6.27) 
Define susceptibility af constant strain Wiss accordiny to 
f $s 
Wi Xin = 4; (6.26) 
q and multiply both scdes of (6.21) by t.%m; to get 
q 6 We Ee = Py t EWN: fiery 5 (6.19) 
ey @ Ny EE > Fn t Fo Mj Tike ke 
Ey Then elimiumate P between (6.29) and (6:1), with the resuit 
ee < : 
. 25 Di= “EN; Fics Mee +6 (6,;+ nije; . (6.30) 
4 Preordingly, Comparing (6-30) wit (6.6), we find 
' s Ws 38 
3 . ee ; . 
4 Cee =~ oj Fikes OF Fae Kaj Pines (6-31) 
3 Jess $ $ _ ,-!,8 z 
E a 6565; + Hii), or He F So Cj > Sey - (6.32) 
4 To tind Coiled in terms of cyt , Furst substitute rhe expresscon for 
Ma 
FE Poy quen in terms of Ey and Si; 1m (0.29) into the sccond cf (6.24): 
4 
a 7; . $ S¢ ? ; 
E bees fnij \ ome Ey - €o mn Fane Sua * Cope Ske ; 
E or, from (¢.31), 
3 
: : P s 
s Ts (cigs ~ €o Nem F oni fare) Ske 7 Cxcj Ee (6. 33) 
. = 5 
“a Comparing (6. 33) with (6.5), we see that ‘ 
“9 € _ oP s a 
ei Ci gee = Cages ~ Fo Yan Fonds fink ea a 
me: or, From (6.31), 
E “ts \ 
Cink = Ciike + &, Kw Cmij Cukl- (635) 


With the ard of the foregoing expressions for the constants ais 


eh Rae es, BE le Seas 


2 eee eee 
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Coke and Fives we may reduce the eq uatiens of equili brian (&.20) ana the 

Comsts tu five eq catrene (6.24) to te classecal form of ceuations cn U? 
Vv 

and ¥, Given wn (6.8) and (6.9), by elimination of ean Fiesty to solve the 


fist of (6.24) tor P, recall that the reciprocal of ae 1s bei; . they, 


upon mul tilicateon, we havt 
a) L 
Pe 7 & i  (F; r Figg a) : (6.36) 


Substitute this for By im the second c& (6-24) te obtutn 


6 ee ae 
een Cetin Ean * Fp Dog) * Cijet Diet - 


From the sccond cgus librium cquation, ES a , again omittng be. Alse, 


Sp: glu, FH \. Hence 


ty bid 


_¢P s sekee tye 
Tape = (Coke 7 oY men Fynij Fait) AKte ~ GoMam fucj Fei 


Then, from the first 2 ¢ (6.31) and from (6,54), 


age 


yi 
which 1s the left hund side of (@.8), as required. 


Next, €rom (6.36), 


wes 


$s t 
Pot ENG; ( EF * fipe Sei,i) 


ae ; 
=~ oi; Hej - Fo Mi; Fake Uh, Le 


Hence, the left hand side of the lastof the equili@riam equuticns (6.20) 


become 5 


$ $ o 
— Eo ic Le OMS Fake Meee” CHa RA : 
ov, from the first of (6.31) and (6,52), 


noe ~-° _ 
~ © +P; a7 Ooke Uk ee egj Cis 


an wae ne nse es ea ke nha AED “Ini AA REID berh AR ADAAR DLN IRE OAS SY OF 


oa: a Ns ie AO et et 


ke qewg 
er ke ie spe cele ny ce 
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7, Polarization Gradient 


Toupm' s form of the classical equations ef clastic diclee tres reveals an 


omission tn the classical theory. His equation 
Ep -4; rel =0 


rg the “eg ua tron of intramolecular force balance’ [4z]} weer dérven ecu 
fundamental considermmtions of the Cull Orie ot electra acces. bat which does 
notavypear in the usual formulations: Granted tHe valrdete of rhe cquateon, 
IT 8 segue ficant that no boundary cond: Fron 1s associated with if , Whereas 
thercis an eg uilibrium cq uation associated wrth eachot rhe varcubles U., 


gy und F, only us and ~ arc accompanied by boundary ends Hons 


fherc 15 no coe thictent of Gr; in the surtace ia tegral in (e.10, te cumyle- 


ment that in the volume vateg ral. This lack can he truced buck to the 

absence ofa Functronal de pendence ot wt on the prlarteaticn oy radient’t « 
B 

Lu tucf, iC we were te start by a suming de penden Be af Wen the dis - ' 


placement aud polart ration and ther gradients and teuncautcia Fler 


the berst gradient Pr would remain along with Se and P2, Calg ay 


and the ankisammetne perl of us would have to Cérdiviavdca -- oa the 


qrounds o f nee etree Frunslationul aud rotakenal tnvuriauce cfwe, 
ihe culy paste fication tor discarding age would be its possiole lack 
of importance Judycd on practical consrdérations. However us witl 


be shown, Pi supplies icems Found ia the (Ong wave limits of Ecmite, 


di ttertuct cquutcons of lathece 


theores of Crystals and ulso Cxtends rhe 
a : 


Coutimuum theory to accommoatée observed Physical phenamcad, 
To éy tend Jouin’s Yara tronul primiple ty account for the contrio- 


ufion of He Polerizatcon gradieat, (tas only necessary fo rc lac€ (b-I2) 
with F 
Hew'( 5,; PEs Ma be Give + ei P.. (7.1) 


Ln addition, we sholl include the kinctee energy So thal (6.13) becomes 


} t 
poem 3 " = 2 TES £ ia BO A IIL AC $= 
Netosta birdy cos (ar raat bE SE a as MO Wd MALE STATENS on abe ax, Se Sar er 


POL RAPE VO WE SO LE PY 
= FS I TORT Te TINTON ETS bist 
gS) ce Se 3 oe oad WEEE ee i Eee OPE 
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‘ae 


ee FOP eI IO EE Tg gh Lite 


HL ASTing eee BS A acces 7 ae Re 7 


Sea MANES DIAN OSEEY POT ERT OEE BF IHR PORT ARE LER TE TD AS REEF VEE A 
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40 
: ‘ 5 , : ! : 
! £ t t e eo t - . ° u ‘ ‘ i i. ‘ ‘ “ 
f at Bedat, wav {ae (F.54ir eR IAV +f tures | = 9, hal 
t, ye ts Vv $ s . 
ter independent variations of Uz y und P: between Freed limits ut tomes anal, 
the only addi tinal creralians As volved urt & hew ae i ailee, _ 
] ‘ | . : 
' “oot (- ‘ owe | ' : 
, ~ Oy See AY) 
i LJ 
} i : 2 a : 1 1° | 

7 , and twa In tearutyons 

: 1 ‘ ; 

: i t, t, t t ; io } 
: Ci edabe| asudesoual | aaeer aad 

‘f 20) uude=s| U;cu;dt-u;du) -| d, bu; de - = G Su;at, : 
ce \ . to : t 2 te /t, . ue 4 

ve oo ’ } ' i 
uJ re, EE ~( ye §P-ds - Ev ds: . 
a Pr 6 ij cv Ns oP) LAC { ncE,, 8P: ds Sybase Et, av. 

with these results and a 16)! tas becomes 

Es : F ; ! ; , ‘ 
8 : f = | 
Bs ‘ = +. r t -° : 
a i at], Liz. t f, ~pd,) bu; +( beri ba E; i te, er. r(-G yt? be]dV i : 

ce 4 . ' 4 

! 


- [ae 66 Se av + [ ("ae) Ut; ~nTe) Su; - ESP, ene(sheyer) belds-0. nay 


a : ' H 
n Then rhe Luler Eg a eens arc ' i 
4 1 1 i H 
BS T, = oo , i] t 
a , Bit GY Fede - 
: , 1 . , v 
q , Cat + &r- ass t ey =0, inv (ho) i 
. 1 ad ' ; 
a 
a oo Hii hye 285 | od 
a : : ; t 
Ri is 
F. , ‘ ! Gi 0 wy V , Ch 6, 
K ; : 
EY ; . : ' : ! , 
S with natural boundary conditions 
Ri . ' n Ty = ti, > es 
i ‘ ' : : 
| ) nvF,. =0, (7.7) 


‘ : parr 
ET Le 


¥ e, i PMP 
MATER se oes 
cb SS Sn into Secchi bse it ae so acres Ripka ete tn aa al, ai ic a tai es 


a(-ebe is Vilees -. | 


Kae 


2 ore 
a + For W*, we teke ; | 

a oa W's be a Ea, P, D Pit 3 Bike He fox t 3 off 0 5s § 
a + Vige Pii Sha * Fisk Pisin . Juin Page : (ne) 
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eee = 


ars x * ey My 
| pA ED Seay a ST ERE OS LR Oe ae ES 


is 
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The superscripts &,P,G, designating fiacd sfrain, polaritation and polarisation grad - 
tent, respectively, will be omitted in the sequel where no ambiguity results . 
Substituting (7.8) 10 (6.19) and (7. 3) we find 


t 
MES = Bi Ph ei ice Pee + Fine Sue » 

. oO i 
Ecj = Jncj Pk + ike Pent ghee Ske * biz, (7.9) 


Tes = fray Pa PSkeey Pak * Cane Ske. 


The field equations (7.9) and (1.6), the constitutive egua tions (7.9) 
and the shuin-displecement relation 5; = 5 lus; eu; 5) comprise te eguatons of the 
augmented Pncersy [44], Its aprarent from the form of the imntegrana of the 
Sur Face integral ia (1.4), whet voundery cond trons other than Ciljare admissible. 
Thus, in place of the traction a, lz, may be substituted the displacemut u; era 
component of cther and the resultant o & the other im the Plane ct right angles, the 
Same possibilities are cpen for n; Ey aud Ki, Finally, either the surface charge 
a (-e by] +P.) or the potential pe may be spee hed, Cf prime importance is the 
fact tout both the potential y and rhe polar; talton , may be spec fied: his 4 
@ latitude not permissible tn the classical theory. 

Cue of the novel proper ties of the augmented theory cf clastec 


diclectncs 1s tts qtcommodation of an electro mechanical interaction even 


for materials with centrosymmetric physical properties, thes ray be scen by 


in spectron of the energy densify W° as given 11.8). For cenfrosymue try , 


SS dine 3 oe cists ate S 


fig? Jijk = 0, since there arc Ho centrosymme tre tensors of odd rank. However, Ars 


still {eaves Soke which is the coctficcenl cf au electromechanscu( coupling perm 


Ene 


re BE by 


in the augmented theery, but docs nel appear we the clossicel theory. 


; 3 : As an example of crestulle graphic cenlrosymmetry, consider the 
i 4 cubic pomt group mm (luternational) 0, (chee flies), the generators fur 
8 which are Ces] 
R 1 0 0 0 9 | 1 0 06 
q Oi Si AO ede On  Pab) e (710) 
ye 
ae 0 oOo -l 0 1 O o -t oO 


6 RGR 
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Applying these to the coefficients im (7. 8) we find 
Tijk = Qo, lis 0, 


a; = Ay Se; 6:57 zh 0Sij5 


POE EAT ORE EE II oe ETE eRe Va Pe Ba 


a bé;. jks + by b:: Out + by (5:4 6 ig * Sig J; a) + bry (Sz, 5 0785 Sie), (7.01) 


RG Pode hoy 


Coie © +12 Sibu + Coal Sig Sie + be Six), 


diinc= A Sige ie Siz bee + deal Sey 6:2 +84 Ste), 


SAPP? RONG TEST 


where 


b =b,, -6,, -26,, 9 C= Cy WC -Lowe 3 d: dy, “di, ~2egy . (7, (2) 


Insertion of (7.41) tn the constitutive equation s (7.9) reduces the latter to 


ee 
ELOY NEE MED TE, CRIS CTA ANG PATON PLES S 


wa eee ee anne ecm ee ene 


y “E: 2a, P; 


3 


€.;? 0 Sine Pek + 62 6; Prt t baal Pi; Pj) ebay (Pee “J ni) h8, (203) : 
A 
t d 6; ski Ske tdi O.. eK + 2dye 5 iyo 


i 

q a : a Seika Poi t fo oy Fok tdgg (P=, = B \+ed: Sine See C2 gy. hx t 244 Sy "6 . 
E The egcerens of motion on ui, P: and p are obtatned by Sabstituring (213) 01 
4 (15) aud employing the stran “Ais placement relation with rhe cesull: 

:: 

i 

q e4; ae ac WA hel t Cyl Were rds; jet Pee Pia he ky dag ( Rai Paul tf ae 

q Gp Urn Saag Hel, rij Pode Paki t biaPa ky rue Pi Pea) (7.14) 
q Tilt 
ry, ° 

E - bar (Bcc Ric! ~ Gy r; ris a =0, 

: hee hee. 

g 
4” In the one chiimeustona! cas €, 

Re 

q : 

a 474, nt), U2euy= 8, PrF (a, th, P-F-0, ¢- ¢lr,t), e3) 
ES 

a | and ia the absence of extunsic Fields f; and cS (Z 4) reduce to 

BS 

¥ > +P, = ou 
fy 2,4, dy 0, P, = CU» 

F a a 

a, u, t by, 0, F “ay P, xO” = =0, (7.46) 


-€ 0% +OPp =). 


a ee a ee ee 
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Ve GEE RNAS eee PORES a MN mT NCA A FR He 


PY iterate ee — enmmment ga mereeenttnmatnne lathe a mein nue neem ne ea Sunt eRe Stee ON nS annem sere pet teneee Ata Semen miter 
43 
To obfaa the disyersion relation for (7-46), let 
(bn -we) an bee 
urhe (3x w ‘ P - Bott we) : Ys Cet Sh wt) (117) 
U pou su shi tu eng (7.17) ta to (2.46) aud clumdatens A, Bande, we offen 
Cub *= pw 35° 
. = 0. (7.48) 
4i§ 6,5" +a, + 6 
qos isa. el Cg tation i hae su that there ere two branches. lose 
deturtencss of the energy density reguires 
au>09, bylu-dye > 9. (2.49) 


Consequently one of the vocts ¥* Is positive real and rhe other 45 Mey ative 
real. Thus, one of the branches of the dispersion rclation ts real and the 


ether ts Ie WG ure - the behavior ut the real branch af low Fregucnc 7 
15 quven by 


lena $7 ]% 2 (1.22) 
goo) Te 
The linet of the lanaginary branch at tévo frequency 1S 


lim & = i/f, 
wb 


(7-21) 
$40 


where ‘i 


2 
ga | Lbutu-Ailéo| as 
(it €gauden 


if the contrcbution ot the polarie ation gradicat ts aed a 
and d,, cere ero aud (7.16) reduces to 


Cu Oty, = 4, > 


a, + op=0, (7.2.3) 
- 6, 0,4 +d°,=0, 
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whick ¢s Toupins form of the Clu ssteal eguateons for the one-disircusconal, (ince, 


ceatrsynmetyic Case. In (7.23), There 6s no electomechanical Coupling ; He 
wave motion 63 dispersronless with frequency given by (7.20) for all g, 
Tf, ugaen , the contribution of the polarizatien gradient ts onmettes bet 


the material is noun centrosymmetrc, (7:23) becomes 


Cu Mu, t ¥,,0, 7 = ea, ’ 
fi ou, + a,P + Oe = 0, (7.24) 


- 62,0 + 2,7 +0, 


so that electromechauccal Coupling ts regained . Lusertiag (161) it (it) and 


eliminating A, bund, we fiad 


gw= Le, - £2 Mayr eis . (2.25) 


Thus, even though nou centiosymmetry induces electro mec hanilal coupling, Phere 
(s still no imagin ary branch and the rcal branchis dispersionless — it the 


absence cf tecon tribution of the polar cation gre dreat, 


‘ 
sad te mS Sli het Selec a Stab at CEs aw 
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8. Monatomic Lattice of Shef[- Model Atoms 


In thes Artile, the difference equations of motion are derived for 4% 
one-dimenxonal, monatomec lattice of the Cochran type lte] based on the Pick- 


Overhauser [47] shell-model of the atom: a core, Comysrising the nucleus ans 


PINE RSE ate ee e 
‘as — ‘6 ey 


Re 


inner electrons, surrounded by a shell of outer clectrons. [he elec.tronce 


es 


polaritation 15 proportional to the relative Aisplaement of the core und shell 


ayo 


of the atom . rTu add, tion ts thes intra-a tuimec rate racteon, acount ts taken 


ot imteratome interactions between cove ead cove, corcand shell aus shell 


} 
aud shell of nearest neigh bor atoms, Ltts shown that the € guutions of 
the lattice have, as their contiauum (omit, the augmented Zyaurinis [48] 
described ta the preceding Article rather than the equateons of the cleswcel 
theory of elustec dtelectecs, The addibronal effects associates with the 
new material constants bike wad Sau stem prumarcla From the shell- shet/ 
interactcon . thes interacteon ts known ho be important A rhe tee tele ime, ef 
lattice dispersion relatrons to dispersion datu from neutron diffraction mus 
urements «af short ware lengths (2c, +8]. 

We consider a Single ude of atous (Fig. 12) extracted from @ three- 
imeustonal monatome laftce of shell-model utoms . The ufoms arc Situatated af 
K,2na where n ts @ positive or negations ta teger and 4 is the lalhce perumeter. 
the displacements of the core und shell ot the atom af Ksnu are A3ig nated 
by u, and S,, re spectively , The force cons font of the intra-atomec core-shell 
taterachon rs de siquated ba d& aad the force constants of the mterahomec 
Core-core, core -shell oncl sheli-shel/ interactous are designated by 2, yf aad & 
resyectively. 

The equation of motion of the nt* waterior atom (core und shell comoncd } 
(s shtumed ly setting the ihertia of the atom equal ty the sam of the forces 


on its core and shell by tHe cores aud shells of chs two nearest nigh atoms: 


Olu, is a) Sea ~Uy) ty (ung Sn) . é (Sie -S,\ 


i é (uy, Uy) -yluy- Sat) a 4 (5, - Uy) z 5(s, Z Sac) = mM u,, ) (a. {) 
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eeeremeeeenne meter «Recent SARA DEERE PELE DISET IAS TOO A CS A “oe 


ee ere 


46 


where mis the mass of the atom. The shell, alone, oF the nv" atom 1s acted / 
upon by its own core, the core and she ef Eauchot the two neta h boring tous 


and, also, the Marwell, electne selt- tell (ahich occupees all space). Meme 


o (U,- Sq) f ney - Sa) + 6 (Sues -S,) 


os 4 (54-Uns) ne S(sy- net) + 8 apa =0, (8.2) 


SO EE eae UA EE ASA PRN EE NCAR EEG ee MOEN TA 


tn which the tnerthx of the shelf fs neglected, q ty the churge of rhe wtom 
uma Ey to the value, afx,-na, of the Mavwetl field . Wow, E, = -LO8S: na 
nd this muy be Expressicd in terms of Uy , the value of the yoy hential ¢ 
ut Kp nu, an Cxpundiny the derwvatwe in an wntinte series of forwurd 


differences [7] * 


E,= 0°¢ = = (1 mw a” ATY, » ( 6.3) 
where 

: A, €,? (ei-€, Va, 

; A = ( ¢ ~2 Gn, r¢,la® 
2 +¥y €n 2 A (8.4) 
By Pn? (kas 3 Prgr +3 Pay -¢,)/a 


e . 


The poleritation oft the ee atom per anit aren ot He Three dimension: 


a ul latte, ts dettued by 
P, = (5, ~ 44) /a?, (ea) 


Then the equations of motion (&lj and (6,2) mug be rearranged fo the tollow - - 


Sai vera atten 


ing forms 
(i rey r bla Al u, + (pr bate” a P= eu, 


(8.6) 
res)atg' Ayu, + Sa’g*a\ ?.- (ardpa'g”P, Og, 0, 


where p- n/a and 4? 15 the second ceutrual difference operutor (diviace by a) 


e 
“ay SScaigascASoh PFS LIT AIN SANGIN Rathod eal ARIAL EPA ADRS IO A OSCR OA A I * 


as defmet wr (2.3). 


AH r a2 2 Pee an) oe LAN a EP ea TA en SPOR AL esa Oo 
POOR 8s, STE PAN EL SRT Sy DY PERS Ee 8 we ee mn Oe ete Rew eee 


Pee EE Sak Oo rare TE LOAN OL C E  OREAE 


T+ will be observed that the differeace equations (8.6) have the same 
form as the first two of (ae), which wre the one- dumensconal as thlereatral 
Cg untion of mottou of the augmented the ory of elastic dielechece . Axcord - 


gly , iT we set 


ae ner A ate RN DS Rng 


\ 


a, = (Ar2yp)a? er 3 b,, = Sak” | 
Mt rey 4% t v4 (@-7) 


Cy 2 (@+2y+6)a”, dy, = (y+d)a‘e” | 
aud emplou the evpantcon (2.1) of the secoud ceutrul difference im serves 
ot derivaliwes, we find that the fust two of (7.6) ave the lowest order 
Continuum Lypod im efron to the equations of motion of the atom aud tts 
shell, resye etively ; 

As may be seen from (8.7) the polarization gracent terms (Those with 
coetficrents b, and d,,) in (1.16) stem from the shell- shell and cuteratomec 
Core- shell interactions — (denti tred by the force constants S and yp, respect 
twvely. La fact the form otf (7.16) «s preserved «€ the interatomec core 
shell cnteractcon «ys omrtted gut the shell- shell iuteractcon 1s refucmed fr 
both ot these mferactions ere omitted, the continuum approrimatcun 
reductes +o (217), whieh ts Tousen! s form of ¢he classica ( eguatcons 
(1a the one-duneusconal, likear case), The second oF (1-17), Toupiad 
“equation of titramolecular force balance", dees net aypear ia the tradi - 
tina ( equations of elastic electrics, bufwe see Hatt 5  Funda- 
mente | equtlioriam condition: the \gui nore om of rhe shell (ce. the 
cuter electrons) under the actin of the core of the sume chou and the 
surroauding Max well field. the COrTES ponding eguatsen in the 
eytended Heory (the second of (U.16)) incledes, be ald tion, rhe ackeou 
of the adjuceyt atoms on the shell (1a tae continuam a pprowmarres) ; 


We have yet to derive the lattice counterpart of the third of 


xs 


3 ales Ladisrsion Eines Nes DiniC ASAE aor n Pamemnaane 7 mM 


(ile) (the “charge equation”) aud estuglish the boundary Conditions. 


Its ‘Uurmacacetrag to yeach these results From considcrations of eucrsy. 
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By analog with (6.12) and (7.8) lefus teke, for rhe electric cuthelpy 


ot the one-dimensional lattice (perunct volume a?) 


HE [bae,4+46 ) Pr + $b, (A,P,)* + ¢ (Osu, 


+ dy (A.PXA,u,) - £61 2"y,)” * Py (a%y,)). (8.3) 


We have to find the derivahves of HH with respect to u,, P, and ,, The first two 


ave Convenhonal and shatghtforward« For Cramp le, 


2 a 
24, 4,4)” = 4. (Umer “Una) = AU ee -Um) dunes - Usa) P 
Un a DUn a dup, 


Z (ner ~4,.\ 5" -6™) ; 
a> 


7 SU Aner Uner Un) = 2G n , (8.9) 


where se is the Kronecker delta . To find the derivative cf O° with respect to 


h,, note first that, from (8.3), 


ore, sa [ asi a 2 CPne2 ~2 Fao t Yn) 


. 5 (eas “3 Crna +3 Caet -¢,)7 seals 


3 Then + 
* od" ¥,, eye mlm mr me m 
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o$(P,, ~3P,.2+3P,., -P,)- nae) > 


2-dP,, (E12) 
where . be ee 
oP = 2 ma A. P, (8.11) 
j ana 
} ALP, = (P.-b a, 


A. A, = ( P, - 2P iy r P,al/a’, (8.12) 
AY P.# ( P. iz 3Py +3. 7 P,.3)/a%, 
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ie. OP, is the expan ston of DP in au infinite series of backuard differences [7]. 


Finally 
+t 2. +, 
; a(d ¥,) 227e, r1¢) Goa) 2-279" oe (8.13) 
oY, De, 


We now tind the € uations of mofton: 
Pr) fase! 2 2 on ef 
- > ix = Cy BU, + dy ALP, = PUn » 
fn 


: oF = Ay Ug + by Ai Pa ~ 6'' Py 2° =0, 6.) 
n 


- ae ~€,0 OTH, +d P, =D. 


The first two of (8.14) are the same as (8.6), with (8.1). This 
justi fles the assumption (8.8) for the electric enthalpy. Also, we have Found 
the third of (2.14), which has the third of (7.16) asits continuum form. 


T# the lattice t's of fiulte theckness spanning an odd number of atoms 


with the end ones at n= +N, the conditions for free boundaries are 


rye 
a ac = Cy AyUsn 6 A AsPey = 9, 
_ _oF = di Alen + bude Pay + be =0, (8.(5) 
D(A, Pe) 

- es, = €,27 Pen = Fog =0, 

d(d* fy) 
where 
Ds fau> Feaen fall, OM yy = (PTH) cane @.16) 


Tu general, admissible beundary conditions are: the spe cification of one 


member of each of the three products, at each end, 
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The dis persion r+lation for (8.14) (s obtained by inserting the 


function s 


Us Act tna-wt) P= Bei bne-wt) 4,2 Cet Hna-ot) (8.17) 


nS 
ke 


tn the eguutiuns aud eliminating A, Band C, with the result 


2 2 -2 
2a C, (I- ¢eska) - gw 2a “d,, Ui-cos ea) 
ut ? a $ o, D ; (8.18) 


STR er on 


2u°?d, (1-cos ta) 2a *b, (!- Cos kau) rates 


zs 


This two- branch dispersion rela tion +s to be Compurc d wif (7.10) —- the 


alispersion relation for the long uuve aypromugtion, For by <cl und real, 
Z2a™* (1- cos ga) —> £7 


So thar (8-18) reduces fo (118) as expected. /hus, the (ong wave be ; 
havior of the real brunch of (8.18) és again given by (1. 20). As €or 
the imaginary branch of (6.18), the wave number at zevo tregwency ‘s 
oe by sinkg = 435 or a= ‘e, where sinh Fy = s , (8.19) 

instead of fa=taf/k. that’s, a/though the low Free uencg cuds of the 
branches for the long wave appronmation (7.18) are gual iFarevely the 
same «5 the low Tr te vanes lint of the branches for the luttice, only 
the real bruaches are quantitatively the same there (Fig. 13). The 
limit of the imaginary branch fer the (ong wave apprommution ts 
somewhat th Error as the wave number at w=0 is not small. Typically, 
a/s might rang é between laud 2; then the Crror in the continuam . 
approxm aticu Sa-cafh as compared ui tt ta/rk, would range from 


47, to (2Yo — on the large side. As wll be seeu, the Imag eiury weve 


number (5 ugsociaked with surface effects wot fuund ia the cla rsecel 
theory audanu errur th the wavenumber affects the amplitude ot these 
effcuts uud thew spatial rate of decay cate the intercer , 
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The range of a/h From 12 stems from calculations, by Askar, 


Lee and Cakmal [49], of the values of the New coustants 6) auddy, 


2 agen neennent enmeeetneeN neem RR gene ane 


: Cand, Couseguentin, f) for several alkald halide crgstals. They found the 
relations of these constants to the coustauts th a taree-dimensionel Cochran 
146} leftee of Vick- Overhauser [473 shell-model atums im au MaCl structure, 
The Correspondence is not complete because the MCl latter cs deatomec where- 
as the continuum ian brtecle 1 5 4 stmyale Cie. not compound) one. Hower 
the long wave, low tong panes limit of the lathce epuatcons hes the 


same form «as the continuum eg uations iu Artele J, 
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9, Ca pucitance of Thin, Dielectric Films 
the Cxample ot te electrical capaci fance of thin, Gitlectric filmes 

Supplees a@ vehicle Lahogie for Compares ons of CK perimental duty with 
the predictions of the Casscaal theory ot elastic deélectres polaria: 
ation g radren | theory and (atfFice theory . 

Consider a metal- drelectnc -metul sandwich whose mradle plane 
Ka-%y, 1s « plane of geometrical and mafertul symmetry, If eud effect are 
melected , the fields are one-dimensional aud the eguatrons 2 Feguclibrinn 
are given by (7T.le) with a, = 0, 

hetus firsl solve the problem wifhin the framework of the clussice/ 
theory. jheu 6, aud ch, are tere aug the tg uatrous of cgice le bres reduce 


to (7.23) wth & = 0: 


a 


ey dae, a, 0 +0, 0, a + OP ~O0, (9. 1) 


Fer a dtelectnre la yer with tructiou- free surfaces af x,~£H, on whch 
are impressed volfuges +*V, the Goundary conditions, according tu the 


classical theory, are 


Toad. = 9, bes cise =4V; (9.4) 


aud the solution of (9.4), subject tos these boundery cou difeons | ts 


/ 


except for additive constonts in u aud Y 5 
u,=0, he yVh , ce Wa lh, (3,3) 


wheve , = 6a, ‘ ts the drelectne Susceptibility. 
Ihe Cee pacitance (per unt arca] ¢s the vato of the surtuce 


charge Cywer unt arca) ty the voltage dro acruss the lager: 


cw ee Fleas 2 Sellen) oe (9.4) 
2y 2h ah 


where € 15 the peemittivify of the chelectric . 
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Thus, according fe the classical theory, there ts no strewn , the 
polarceution 65 ur Korm , the potential varies lineacly through the Hecok ness 
and the capscctance Is inversely proportional to the rhickness , 70 Rata 
qraph wf wnverse Capac! tance vs. thickness ts Ze straight ling through the 
origin, However, ina sertes of Ekperiments witha variety er vere the 
dielectree films between metal electrides, Mead [59,5 : s2J Faia bs 
dif€crent relation betwoen (averse Capacitance aud racchuess, tee ek: 
perimental data Tull on stracght liaes whieG if exttaded fo eecre 
thickness, have positive intercepts of inverse Camuuitadce, us illastat 
éd in Figs t4. Tucteatly, [sol], Mead Sugg ested Hat the ancnle ly night 
bc due to penctral ron of the field tuto the elechrudes: but supegedntly 
{513 he ehandmed that view, although , meanwhile, (T hu Yee sup- 
ported by Ku and Ullman [s2], Au alternate € xylanation (s found . 
in the augmented Heorg of elashe drelectnis [54], . . 

The soluttan of the augmented equations of egucli orice , (hte) 
with u,> 9, analogous tr the one just given for the ia ae Cquattons , 
requires a boundary condition in addition ta the surfuce traction and 


potential, The additional conditrow con be the Speci ficatron of He pote 


ization atthe surface of the delectne . Vow, the polari ration ata 
boun dary « ve the dielectnce twa mefel- dielectric ~ me tee Sand wich wil{ 
depend on the phy seca / properties of the adjacent efectrude' aud 
metul- dielectne tuterfaces aud these properties ure cutside the 
Compa 95 ot the theory of dielectrics, towe ves, Stick the pelanh cation: 
in the meful ¢s zero, if is veasonuble 1 assume that the sae face paleriee- 
tion in the dielectnc, 1f rot actually gery will We between wero and 
the classical value given by the second of (39.3), hus, us ming Hat: 
the twu electrades and interfaces arc the Same, ther influcue on the 


Surface polarieution may be infreduced phenomenolug ically , by >eThing 


the boundary Couds Keon 


Tee wg ee 


ss. IS 


ea wo ale 8 nates te 


pea deeee ng we 


[Placer TT KEHV/A, = DEKE, (9.94) 
The: 'elassteal Ssnatbiss 5 k= | while k>0 describes continuity of 


polariautio n actvoss the inter faces. ‘We supoose also, us ts wssumed iy 


the classical the or that the traction weruss the wi rere \S e€Vo, 
% 4 


Evom C415), this condition es 


i 
: ” ! 


Rta lei = La ane ai e\as 7. 858) 


Finally, we suppose thut the voltuges applied te the. dtelectrie at the 
cn tev tu ces are again 


Lelaeh - £¥. Gee) 


s 


We have now to solve (7.16); with ci, = 0, Gabjcet te Phe 
Youn dary conditions (CY. Sa bc) her 
ee B, cosh * 2 P= Az +B, cosh Me, wshysx, +B, sinh * ; (9.0) 


Upon substituting (9.G) in (r16) ve Liane. ' 


* ALEMAS equals = Ley Bs Leads | coe DY 4 


; Ee 12 12 ) 
yi = { buty-¢ Bit tay ) - - (aes Fy (by Cu - 4 a, \ 
\ Clay rey! } c(t 4" ) ai 


where 


i 


te. £ ts the SUME a5 in Ci22); 


4 


As for boundary conditions , (9.5 6) ts satys Fréd i deute- 


cully white (9.5a), ana (9.5¢) become . 


, | A, + 8, cosh (h/£) = -ke,yV/h, 


. +f 4 (9.9) 
Ayh +B, tah (h/() =V, 
vespectively. Fvau (9,7) and (9.9), 
Ay2 (8/2) cosh (n/L) + kh, : 
CY. 40) 


By = (I-k)yV/ [4 sinh (h/2) + Ch/E) cosh (W/(), 


peta 
2 ee = eee een: od 


IF 


The remaining constants are obtamed from (9.10) and (9.7). 


The capecitauce ‘$ 
le oY” een € | +(ky£/h) tanh (h /e) 


2V : Zh | + (4 £/h) tanh (h/) wy 


C = 


rF we igneve any voltuge drop that may occur in the electrodes. 

Ta Figet©, the curve marked k= 04 shows the relation be- 
tween normalized inverse capacitance aud norwmaltzed thickness, according 
to (9.1), for the case 4210, K= 01, Gleulations from Meads date, 
For small k, indicate that the materia/ property 4 for all bis dielectnes 


ts of the order of magnitude of a few angstroms ; 


and thes ts sap- 
ported by Askar, Lee and Cakmukes calculations for alkali’ halides 44], 
Hence, Meed’s data, which do not extend below a thickness of 30 A, would 
be well to the right of the knee of the curve k =Ol,in Fig. IS, and so 
would give the appearance ofa linear relation whech t€ extended yu zEerc 
thitkness, would have a non-sero postive, Mtercept of taverse capacitance. 


Thes intercept, according to (9.11), ts 


Cy = 2h (I-k)a/e, (9,42) 


Te kel, the mtercept reduces to ¢ero aud, in tact, the whole solution 
reduces ty the classical oue . However, it seems unhkely that te presence 
of the metal would not influence the polarization of the declectric at the 
metul- dielectric interface. 

The polarization and potenti / vary across the thickness of the 
dielectric as shown ih Fig. 16. The absolute value of the polarizateun 
is almost uniform across the major portion of the thickncss aud 
slightly less than the uniform polarization of the class (cal theory ; but 
then dy ops sharply, Near the surfaces, te boundary values K times the 


Cla ssical pelarization, as syece fred « The poteatia/ hes an almost unttornr 
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gradient, less than the uniform gradient of the classical theory, over 


most of the thickness, bul then increases sharply on upoproaching the bound- 


arce 9, The extremely lecaliecd surface effects are accounted for by the 


coutribuhou of the poler(aateoa gradient the stoved energy and «are 


associated mathe mutically, with the imaginary branch ofthe disperston 


relation —- a branch not presentin the classical Meory, 


Vew consider the analogous solution of the lattice Cguution s (8.19) 


with aye 0, subject to boundura Cinditians Gualogous to (9.5)+ 


wenheny¥/h, CyAay td ArPen =0, ye eV, 
where h Na. We tuke 
Un? Bicosh4d 5 P= A,r Bicosh™* , 2 Ayna +B, sinh te 


and note that dF (na)]= [2,F(«,)] and 


Xena 


A’ cosh 2 an =a*[ cosh tHe teo hor aay nee sh A. = Oy “sinh th cosh ht . 


Then, substituting (9.14) in (B.l4) we me a Waleg cusses CFU, 
A HAS,  Bl= ABs /e, =~Acu B /eodu 


aud, in place of (9,8), 
sinh (@/2A) = a/2k, 


(9,13) 


(9. (4) 


(9. ¢5) 


(9. 16) 


(9,17) 


where £ is the same as in (9.8). Thus, C9.te) ts the same as (9.7) with £ 


reylaced by Ay but £ and A ave not guite the same : being related im 


accordance with (9.17) aud differing by some 4 to 12 percent, ws remarked 


following (8.19), 


Applicaton of the boundery conditins (913) leads, by the some pro 


cedure as for the continuum, to 


o_ Ci-k)4V ‘ 
? 37 4 sinh B re cosh a 


$5 at ae oho Dons toe 
RID ETRE RAE EO We 


(9.48) 
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s7 
analogous to (9.10)+ Finally , the capacitance 1s 


€, Oey ~ Pers . € be CkyASh) tanh(h/a) ; laa 
av Zh (yA/hn) tanh (b/d) 


Cc = 


Thus, the euntive so lutcon is identical with the continuum one excepT 
that the displacement andl polarization have sig ai ticoned only af 
the atom sites and 4 Is replaced by A. Tu particular, the curve 
4 Fig. IS 5 applica le to the lattice if £5 replaced by Ath both 


abscissa and ordhate. 
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10. Surface Energy of Deformation ana Polarization 


4 y Tt will be observed +hat the in terna/ enersy deusity af de furm- 
¥ ation aud polarivation , (7-8), Contuims the term b: P. « dinedrein 4s 
a polart zation gradient. As remarked in Article 3, the analogous tarei 
q linear iin the strain, say 5a; (1S omitted as itis of no consequence. 
4 It leads to a homogeneous stress which in « bounded body with a free ; 
Ee surface can be removed by a homogenesus strain which tn tern, can 
7 be vegavded as the reference Contigurateon, The situation 1's otter - 
q wise for the term linear in the polarization gradient, Al thous this, 
q Foo, produces a homogeneous field, itis « Field of Ey; rather thay 
E Te) , &5 may be seen by referring to the coustitubve equations (7. 9), 
q The removal of nN Ey; to free a sur fac: (see (7.7) results, ass iilus- 
4 


to the cutenor of ~the hidy, 


I 
4 trated below, tn a polartzateon and stran which decay exponentials 
As is described ta Article 3, the energy regucved to sevarate a 

body take two pacts Comyn ses a boud energy aud a remainder. The 

former ts what Would be required to break the atouc bonds across the 

surface if the straw, anduow the Po larizatcon were prevented trom devel- 

oping, say by Aypothetical mechoncal anal electrical exterual tields. the re- 

3 lease of the hy pothete cal fields would rese/f ia ce deformartton quel polariz- 


atin, localized nearthe surface, with which ts associated a su rface energy 


of defor mation and polarczatian — always negatwe . thes, the euergy reguwth 


a 3 2 
§ +o separate a body tate two perts is the boud eusray less The ebsolute 
aim c 
y Value of the Surface euevay of deformation aed Pelanisefeo,. The later 
4 1s introduced by the lineer term, 6; ae in W* aud a formula For 
‘ that part of the surface ene fgg nae be found as follows 


The Tolal energy 1 the syste, a eguili brent, ‘s 


[,wav= [Lie +seg. gldi = [wiv + ee 2% E. Fi av, (10,1) 
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where V ts the volume occupied buy the diele chee solid bounded by @ surface S 
separating Vi from an outer vacuum Vis and vtzyVrv' 


By means of the coustitutwe equations (7.9), the form of “W* may be 


converted from (7.8) +o 


be T:.§.. pis ~EtP. +b. P- 
awre Ty Sj rE, Re E. P. + bj re (to.2) 


Through the use of the chaw rule of differeutiqttou aud the divergence Heovem, 


[rz sya = [1y ae d¥ = (Ursa): Ty 4s] av 


jst 


= [mG 4 ds eB Gi a, 
Esse = [LE h)a “Ey, Gla W> mj Pi ds fii BV, 
ois n. b;,P: dS, 


[ite 4 Ver [(y; ¢),; °F: eui¢ldV 
fn Ce; i 4ds - ie Cu dV “f.€ 9, 


where Cyd ts the jumy iP) Y 


-qaeross S, pou lin ser ting these vesults i 
~ lot aw yt 
(1v,t), we have 


fewedv= gf nileP+ Tea rb P+ eledylds 


ee, - ; 
aL a apr Eh + EOP rege) dv - [eeupav. (23) 


Application of the € guili rium equations cee. (bi )and (7.6) with u 0, 
to 3) reduces it to 


[we = Ef (tye; r6°2;) )av es ( nel Tayu, +6; 0; Cesk ea] 7Q) +6,2 JS (10.4) 


ev yd 


Then, a the absence of external fre lds (fF; 2E*s 0) aud with the boundary tree, ve, 
(47) with T20 (10.4) redues to 


(10.3°) 
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This is the energy of deformation and polarization which must be added ts the 
Vond energy te obtain the total energy regucred to separate the materce! 


along the surface S, bs sndicated by Schwarte (sl, the additonal 


eherg ty always nega tive a5 a couse guence of the positive de Frinteness 


of the quadratic part of the Energa cleust ty which, tin ture, 15 reguired fer 


STIL 


stability — t.t. because energy must be stored rather than generated, during 


the applicatren of external body ov sartace forces. Thes, 
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wv ={ls 66. + (Wi=by P.\ldV + [bi Tele 


i [34; ¥ +(W'-6, PAV +f nzb:,P 5 f (0.6) 
Ee Suvtracting (10.5) from (12.6), we have 
: [Laqeei(wt-e er ddav #3 cb; ds <0. et) 
a . . . 
4 But the volume mtegral m (0.8) 7's posttire . Heuce 
a (nz be Phas < 0. (0.8) 
ig S 
Finally, we may define 


bese (10.9) 


as the sur fuce cuergy of deformafron aud polarization per vuit area, 
sometines called the surface tenscon . 


There follows the solutiin for the dioplacement, polarization ; poteata/ 


: anc: Surface energy of deformutton and polarization 1 the half - space of 

7 a centrosymmetrrc cubic crystal bounded by a free (too) face. for Hus 
: provlem the fields ure one-dinéusional and in eq uililoreum 20 the! for 

: the half-space 420, the equatuns ot motion (Kis) veduce to 

udu, + dys = 0, 

q dy D.u, + by WP - ay ~ BY =0, (12.00) 

4 — ete + 9,P =0; 

: 
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and the boundary condi tiwas (71.7) bewme, on X= 0, 


C1 %U + dy OF, = 8, 


: d,, 4, + bi? 6 = -b., (10.1) 
- G0 + i a 0, 
Consider 
u, 2A,e x, /2 ; P, = A,e , Y= ayer, (10.12) 


Upon substituting (10.12) into (10.70), we find 
A,=-A#Ale, = Leu Aleodn (10.03) 


where 4 ts ugain given by (17.22), With (10.13), the first aud third of the 
bounduvy conditrons (lout) ave satis fied identically and the second 


boun Aury condition yields 


° A ei ks bedu 


= on (to.4, 


Cy A (au : 65") 


Theu, from (10.13), 


bo by ers) 
A, = ——2— Poo ek ere. ONS 
4(ay+ 6") ° It ay €y 


I+ may be seen that the freeing of the boundary results in a freld 
; of displacement, polarization and poteahkal lacalleed at the surtaces dewuy ing 
4 : tuto the interior with decuy cousteut Lo Associated with the localized Freld 
's a sarvtace energy ot deformution aud polarization, Vv sack area, 


° given by (o.9): 


Ah ARISE irentnatannccereamimmnarene =o" 


ot 
ws=- . ‘ (10.16) 
AL ay (it am €-') 4 
4 Note Hat, from (6.26), (6.28) aud (6.32), 
a: rae 
; ; ae Seek ey (10. 17) 
Ef Ce. the dielectric constant in the Lido] direction. 
a 
3 7 SXaue 4» “<P A 
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As noted at the ead of Article 8, Asker, Lee and Cakmek [49] have 
made calculations of the new material constants for several alkal: haiide 


Crystals. From these they find, i the case of Mall, for example, 


fe 1.3 x407° cm, WS. -59 erg/em. 


(erg) 
The decay constant £ 1s about half the distance between nearest neighbor 
atoms in the Macl crystal, indicating extremely rapid alecay from the 
surface into the interior, The surface werqy of deformation and polari zateon 
is ahout 30% of the bond energy , 50 teat it is Car from negligible. Bork of 
these vesults conform tu exwerimental observations [202i]. 

Additional solutions for surfuce cuergies of deformation and palers- 
zation (in fsotropic matertals) have been obtaned by Schuwarte [54] aud by 
Askar Lee aud Cakuwuk [55] for infernal spherica/ aud cylindrical Surfaces 
and a crack. The latter have used average values of ther previous Caltu- 
lations of constants for alkali halides [49] to obtaca some guanteleteve re~ 
sults. They td that rhe absolute values of Ware reduced by curreterc of 


internal sucfuces Gud they algo corclute thatW'es ficte at the vout ofc 


creck. 
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Ul. Acoustical aud Optical Activity in Alpha Quartz 


Acoustical activity (rotation of the direchon of mechoarcal 013 placement 


alony the path of a transverse, elastic wave) has recentla been observed by Rue 


[56] m X- Quarte . The possibility of rhe phenomensn appears firsk tohere been 


=m ine Es 
SREY EET GETS ITY UAT PR CARD LITT FPG SF 
= 


mentioned by Silin [ST], if wes accounted for by Toupin L19} on the basis of He 
first struin gradient theory (Article Sit). Partgal and Burstern [58] found ou egui- 
valeat result by assigning dependence of the elastce stiffness ou the wave vechr. 
Lu the presont article :t is shown IF 9] that both acoustica/and aptecal actuitis 
are accounted for ta the polarttatran gradient theory ( Article 7]. du what 
follows, the field cquations are exhibited for the coupled clasts - clectric - mag - 
netic system, the jrroblem of shear waves along the trigonal axis of K-querte 
. is solved, formules are ob facned for the optical and acoustical rotatory 
powtrs and numerical valucs of the new matercal constants a the formulas, 
are calculated from experimental data . Essen fealla, the theory has (t that the 
appearance of opt cal activity depends on an interachen between te olariz- 
ation aud she poler(e atign gradient ; the appearance of acoustrcal activity 
Geyeuds on interactions of tHe stracu wth both the polarisation axd the 
palarization gradient, and is aysent if ether uterachion is missing : 
wheu +e. magnetic field is taken wut account, the Gpuatous of 
Article 7 become 


Tye t= P4;» Cit) 
E+E; + E+ EF = 0, (11.2) 
' bik + B; = 0, (11,3) : 
Ai Bin aj SB ~ P=, as : 
f é, Es + Pi =, aur) 
&. 20, Wi.6) 


? 


where B; 1s the mag netre F hax deusi ty, peo is the magnetic permeadilvty of 


a vacuum and 6-5 ig the unct alternating tensor. Also as before, 


Bibeapuistecd esp baiideare st Aes 


a 
= 
2 
xt 
® 


NLT REESE ESE ETL REIS APE EE Ae rte SH GOL ENE BE RENTS Are en AE Mae Sa 


6+ 
t t t 
1,22, eee 2We F Ey, 2 aw (1.7) 
1 O8,, dP; aR 


where 


Ww 5 aes PP; P. + sbeiKehie Pere #5 Cts ne Dey See + iin % ee * Fain P. iSite Die Pe Ph; 5 (1-8) 
Sar (Yi rus;). (i.9) 


Ir ts convencent to eliminate 8; 


curl of (ite3) trom the time derivative of (4) with the result: 


at the start — - by subtracting the 


Ee EL iy = MoE tee? 3. (14.10) 


Equations (iS jJand (hed ave not independent cf (ut) and (U3), respectively, and 


may bc disregarded for the present purpose. Thus, (0), along wet (tt) aud (td): 


(ut) 


™| 
+ 
th 
+ 
aoa 
+ 
™ 
“ 
~' 


an = Evi = €, Ho E. ta P 
are the fielé equations governing mechentcal and electromagnetic waves, 


coupled through the constitutive cguations ! 


: Pp ; ») 
Tes = Coppa Se * freasPe + Facey Phe » 
2. pk 
i 


“ 


Fike See + kk t jjue Yee (tht2) 


ey = deine See Uke F, whoring Pek Pw) 


which ave obtuned from (U7) and (tg). 
We consider plane, transverse waves propagating along the Xg-axis | Ue. 
Uy, %; and Cy are tero and the remaining 4; , P., and E; are functions of x3 


andt only. Then the field equations (41.1), with f, and E° zery, reduce Fo 


Ti33 > pu, Tes, = Pas 


+r &, 20, Ey, E, +E, 79, (41.13) 


Ey y= €& Ao E, My F, ’ 


Ey 


Fass 7 FoMobat Moles 
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a era we ee Tete eo 


FRE SL EN CT RT EEN OEE EHEC 


an EES RTOS AGF UAVS SIN oO oS ETERS EN I PO THI ER RI RE TEESE" 
me : 


Tre einer es or NR ee 


65 
and the constitutive equations (11.12) reduce to £593 


oe 2 _ Pe 
Ty = Sea Ms ~ Fa Pe tre Pi ys Ty2= Cae Ua,3 t fia? + tie Pys 
~E, = fig U2 66'X,P + ju P3 Py -E7=-f, Ast Getih-jnhs, UL 14) 


Ea = Sag Ui — dir Py Oss Ps Esr= dag latin ® tose hs 


for the crystal class 32 (Internationa!) or Dz (Schoenflies) [60] to which «- 
y 3 
guerte belongs. Ln (ete), the abridged notation 


a HL, 22-72, 33-73, 23-94, 327°7,3817°S,—" 138, 1276, 2199 


1s used for the subscripts attached to the material constants . 
Luserting (4.04) iw (a. 13), we have 


P ° oo 
“oa 4133 ~ Tig Pas td7q4 P33 = pus 
P e oe 
Coe Yass * Tia Ps tre P33 * Mas 
ny -t . a 
. D1 M139 72h P23 + bos Piss fig 42,37 & Xu P,+€, = 0, 


(W415) 
41g U2 33 T2jir Ay + bss P35 * Fe 4s “Ku Py + Ez = 0 


? 
E133 = €, MoE tMoP, ; 
Exsy > & Mba ter 
E Vow take 
34 uz A, smg(ay-vt), U4, 2A, cos o(x,-vt), 
y F P,= ©, sing(x,-vt), P, = B cos ¢(x,-ve), (4.16) 


: E, < C, sing (X5-vt), E, =C, Cos OXx,-vt) 


am and substitute these fuuctcous in (ILtS) to find 


(eee ov IGA, ~ Fibs Hd296 By = 0, 
(Ce -pv'é Ag ~ Fig 8, & dag $ 82> 9% 
dag 6°A, ~ £4642 + (oo ee, XB, = 2576 6,75) =9, 
dug tA, fig As t (bp S 7 A, By ~ 2 jn B - Cae 8, 
Mev?G, + (E, Mov -1)C, =9, 
by lo v-1)C, 20. 


(44.17) 
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Samana Ea 1 


oe 


pice fat Ne NOLL tt 


Adding and subtracting these equations th pairs, we have 


ce 


Gee pv" MA, tA,)+ (44.6 + Fi4)(8, * g,)=0, 
(46 £,)(4,2A,)+ (o.. €*+e3'X,, = 2). 6)(8, £B,) -(c,e¢,)=0, (41.18) 
Mov" (8, + B,) + ( €, Mo V7~ IXe,+¢,)=0. 


ss 


VP REE a2 
2 ieee am marae 


5S 


nd ae 


Thus, there are to soletions, Each covves pondin, te crreularly poleriacd waves 


WEES 


vn 2 
oJ * Rare oe 'f Seep rts oe 
yt . ‘t . 
° 


Lol, ¢.222] vince the amplitudes mush Satisty 
APtag Betty, Cpr ees i. 19) 


with ether alf upper Signs er all Icwer signs. Upov substituting (1149) 


tafo (IG), we see that the upper anol lower signs give right und left 


a civeular polarization, respectively . The cquation for the velucities > cb- 
q tamed by Setting the determinant of the coctticients of the wmpali utes in 
F (i. 18) equal to zero t 
' 
p 5: : 
Sag ~ PV din = fy e 
4 bbt fu bsg SE Xy Figg) fe (ag 
0 | we ve Eo) 
| 
4 Ths is @ Quadratec equation tn vy so that there ave two fairs ef aoposite ly 
3 circularly polarized waves. Each paw of such waves Combines te produce 
d a linearly polarised wave with a rofuting direction of polurization Lot, 
¢ p- 222). Thus, we have two cases of rotary polari ection, (hese rae o€ sden- 
BE: tified as optical and acoustical by separuting out first the clectromagnet : 
q Ie part of (W-20) and then the clectra mechamcal perl. [he two mea , 04 


fuct, be considered separate ly owing to the large vaho ef Frey weneses : 
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Cot the order of 10°) at which the two effects arc ebserved 


Bases 


Bal a 


Ei 


Kelis 


the electromaguehic part of the determinant 1m (u-20) 1s rhe 


ea 


minor of the upper lett element. Thus the par of optical velocihes 


ts given by 


Tal 


tte ee ete et ee. ad we EOP cE 2D 


BS et ETL OO RO SOTO Wale of AHN Pao ae SS REESE EE cet GABE + 


de ed RNS) le a DOE Fea ec ee Boe ee ee eh See aoe 
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2 4 . : 
bee St EX, 5 ej 
655 + & Xu * Sir _ =0, Uae) | 
| ms -€ i 
which yields the dispersion formula (él, p.426] . | 
f 
ng b= (xy F 26, jn + bobs 6)", | (wea) i. 
where n, are the indexes of péteachien : ! 
| : 
ne = C/V, (11-23) 
i 
and cus the velocty of light tn vacuo + : 
c(t, us). on >.> 
| 
; ] i 
; From (1.22), we have 
: (ne 1 = (ni =) 4590. oF oS of (2s) | | 
Now, defuse n= s(n, +n.) and assume : mo : , 
[n,-nl« neni. (26) 
; { i 
Then (1.25) becomes, approximately, ; 
Co 
: nnn, =2ln=1)e jh, - 27) 
q . : { 
: where 6, = @/n, 1s the wave number in vacuo. i 
The optical rotatory power in radians yer unit length, Us given | 
by (6, p.222) : : ee geil 
4 
bon = 5L,(n_-n,). . (wgs)’ Gg 
° . ‘ = 
Accordingly, from (t.27) aud (1.28), 4 
2 
2 2 z 
Gop = - (w'- i) €o Jar be (11.29) 4 


: 

. cy ‘ ysl 
Broan ely 
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is the formula for the op tices ro fa tory power in terms of the average 
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index of refraction, n, along the optic axis ' the wave janis tn Vacud, As, 
= ame, t the Derenenrel sn bane & and the edael eouehei i797 = yises 
which , as may be seen in OH @), measures the inferaction he Fvee the pol 


‘ 


arization and the wolarzation gradicnt. 


1 . 1 


The electromechanical part of the dete'rminant in (11.20) ts the min- ' 


or of the flower: right elemen i so that we have : 7 
H 


id 2 ! 
Cag PV dig os Fie ‘ ' 
‘ ee eae |. 0. (1.30) 
buf + £, bes ¢ *€, Xu + 2 ing ig 


for the equation determining the veluccties of the tw acousticul wuves, “Si 
in fluenced by the quasi-static polarization aad Polari eatcon gradient, Frou 


(u.30), 
| : PVs = gm bag F fue) (bag dy Fin). (RL, 


In view of, (11.19). and the inequality of yYandv., the superposition of the | 


two waves yesults tn rotary polarization (acoustical activity) with acoustreg/ 


rotutory power 


ee aan | ae (32) 


where w is the circular frequency Both waves are dis persive. At the | 


: 2evo Frequency (long wave,) Lumet, érom Ul. BI) and (6.34) wth (6. 28) 


1 : ! 


' Hie pve = chy - eefia/ Xu =e hag ’ : (11.33) 
4 


, whieh ts #le result (without acoustica/ activity, since %=V_) that would %e ob - 
tamed ut the contri bution of the poleri eaten Gradient were ome Pte , ce. (F 1 
dies be. and jyy were assumed to be zero. As the fre quence increases a 
zero, the absolute velocity def ference , yu, at first increuses; 50 that: 
the acoustical activity aypears ane tacreases.: With further Imerease of 


frequency, the velocity ar f fereuce again approaches zero, since 
: j f : } 
' 2 . i 


z ‘ P 2 
thie evs = Cte ~ cha/bss » 


(11.34) 
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So that the acoustical activity diminishes ; but His is undoubtedly beyoud the 
range of applica bility of the continuum theory. Up to moderately large 


wave numbers, CUl.3i), with (6.28)and (6.34), ques, to the First order in @, 


ve /Vs = = (dy ~ jar C14) Org Cline > (11.35) 


where vi = Cia /P: In this range, the Frequency 18 Approximately jorportranas to 
the neve number: w=y,'s so that, from (U32) aad (3s), 


%.= (dre “Iq Cra )eie? w'/ (S4)*. (11.36) 


Thus, at frequencies up to, say, 10'°cps, the acoustica lrofatory power ts 
approximately proportional to the Square of the Frequency aud depends on the 


constants P, Cy, and ce which are commonly encountered in preaoelectricit y 

theory, and also on the constants dig and ji7, which control the coupling 

ot the polarization gradient with the stram and pelarizateon, respectively. 

For k-Quarte, all the quantities in the Formula (1-29) for optcal 

rotatory powcr are known except jiz. Thus, for left-handed quartz aud the 

sodium D (ine, 

Oop = -379 vadian /meter let, y.tatl, 

A,= 5893 x10'° meter fol, ptatl, 

nz 1.5539 {e1, ptail, 

6,= 8.854 « 10"'* farad/meter [0], 

ae Jig Gop X; /ate(ri-\)* = 0.19 meter’/favad | (11.37) 

With the value of jyz known, abl guantites in the formula (1-36) for 
Qcoustical rotatory power are known except dyq. In partrular, Pine [56] 
finds thet the acoustical and optical actiuties have op pus te Sigus and the 
acoustical rotatory power along the trigonal aus 3 about 220 radeeus/meter 


at one gigaherta. Thus, tor lefthauded K-guarte, 


ea 


t 
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Ge = 220 radian/meter isk 
w = 210? radian/secone Use], 


P= 2.65 % 10° kilogram /meter® [ozs, 


it} 


Cay = S794 X10? rewton/meter® [63], 


\ 


Cig = ~ 9.0406 coulomb/meter* [63] : 


Heace, hg é 

dia™ Ke (<4) Te% PY + jG = ~ 174 -0.0077 volt. (1.3 8) 
The second term, diy Cage MS negltqeble 77) 6 om parison with the First su Hatwe 
muy drop the deyendence Of Acoustica! rotatory power On Ji7 ot, on the inter ~ 


uction between polarization ane polerisation gradient, aud reyluce (i. 30) wth 


De = dq ei4 POH (Gia) (uh 39) 


Thes, according to tes theory | rhe Presence of acuashica/ actinty 14 
A~Quari> depends on the existence of the pi ccclectroe Stréss constant 
Cie, = Cray, ated the interaction constant 444, =dy223, Sefmcen stram and 
polarization ; Whereas the f réséuce (.) f optical actwity depends only on the cx: 
joteuce of phe suteraction constant 11, =Jisz, be fucen polarizuteon una 


v? lucy cation a teat. 
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12. Diatomic, Elastic, Dielectric Continuum aud Lattice 


i, Poten tial and Kinetic Energies. 

A continuum theory of a ery stul with two polarisable atoms per cell 
may be produced by combining the procedures employed in Articles + und 7 [eal, 
Atoms in the two component continua are again idents Fred by land 2 with dis - 
placements now desiqnated by aoe Lr, and Polar ctons by po as tas - 


trated 1 Ftg, Ue while polarizatcen gradieuls, Strams , pelative displacement 


aud relative rotution are denoted by 


(i) (C7) (Cy) 4 (2) @ # + & 
fa » Fj oes ts rus) gM ee he 6 EES (Mi -42,) 5 (12.1) 


re spectively. The potcata/ energy densita of deformation aud polarization ts 


assumed to be a quadratic function of Ses pel pot “if ce wf : 
‘a 


wt - t ay (a) ) QW 9 a) * + 
Ww WiC Sep Sep y Pein Pe a Pies Pie, Mes wey), 


(12.2) 


The total potential cuergy density, W,ts again the energy density of 


deformation and polaricateon augmented by the cneryy density of the Maxwell, 
electric self- field E; P- py t 


~ “6 
W: Vt 36 Ete. (12.3) 
The associated kinetee energy density (s, again , 
Oe, Ge, te) ‘ 
T#Z Gf ue 4 4, KL25 6=b2,3, (12.4) 


where pe” aud a are the mass deusities of the two continua 


S 


and the 


‘over~POot” notation designates dif ferentiation with respect te Home, 


ii. Field Equations and Boundary Conditsous. 
The field eq uation s and boundary conditions, corrcs ponding TO 
the potenta/ aud einetce energy densities (12.3) and C124), may be 


derived, as in Article 7, by means of an extension of a linear versiosi 
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of Toupin's C42] variationdl principle for the classicol theory of elvstee 
dielectrics, The extension accounts far the contributions of te polerisc Foon 
gradient the two continua and the kinetic energy. 

First, we define au cfectric enthalpy deusity — the energy densit 4 
dim it ished by the product of the Marwell self-Feld, Ec, and the clectre 
cis placement D; : 

H=W-E-D,, (12.5) | 
wher | 
(17.6) . 
: eae (12.7) 


P; =P fe aa Gel, « (2.8) 


The tacluson of the term Gy rim the polarization density ty Was suggested 
: vy P.C.%. lee to accummedate apylicateons to Conte crystals . The product 
E of thé charge density By and rhe relative cisplacement u® yepresents the 
nouic, or atom, polarization «> distinguished from the clectronce polari- 
Zutiwns Pees Tn what follows, + 9, 's not set equ! to Bre He seperscriff 
i=l identifies the positive fon and t~e2 rhe negative (on. 


Inserting (12.3) and (12.6)~ (12,8) tm (12.8) we find 
q «) 
a H= W'- 56,4; fot BilP i +P + e,ui). (12.9) 


Ina body occupying a volume V bounded by a surfuce 5 separ 
ating V from a Valuum Vv, the extension of Toupin’s variatona/ principle 


tukes the form 
("af at{ (T-H)dV+ af “ttf (fdas ESP Erg, Sut )dV 
q +E ode : a A520, w=l,2, Cidt0} 


4 for independent variatcons of u ie Ag aude befucen fixed units ak times to 


and t;. In (z.fo), V*= 2¥+v", aes aud t"! ore fhe external ody forces 


a =e ee Se 
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ang surface tractons on the tuo continua and E, 13 the external electric Field , Mou, 


c. &) Oh ble na W) ey &l te p tet 
LAT PS) EP SPM EL SRM) + Th bul + Tf Sens , (2. 
where 
(e) ywt ® ww ce) We #_ ww *_ owt 
Ly = j@) 9 . = ~ Spt ’ E;; naa 3 in Ts . ue = cane” . (12,12) 
qj ¢ ue ¢ a | 


Mote thet the EX are the eftective local frelds of the two Conta : @s 
distinqueshed From the Mawell self - Field E., 


; 
3 
: 
. 
3 


weth the chain rule of dit ferentiation, (2.11) becomes 
« + t} © * 
Bw'e ELT ge AOI TP Ddu ~Z (Ee ee aR” 
eZ (TS CO HE baz}; Pa z (Es 5P™); 
Treating the remaritong part of H tm (12.9) similarly, we have, nv: 


§4=-Z( 7, ~e) 60 


Uy 


Paci(ty -7/"- 9,8; )] du" - Ze; aad 
+2 Le, Oi ~ Pee - (-1)* 944 ue; 15¢ rE Te CMTS [&u"} 
i E (ens Pe): “Elf ee; -P@-¢ 9,45" dS: $ (12.43) 


and, in V'! 


aM = 6:8 ~ 4 Ce5¢), cn 
Further, from (12.4), 
g(r ee z iy? so Ch) buat : (12.15) i 


Lnscrting (12.13) - (2.15) in (12.10) and ayplying the divergence theorem where 
appropriate, we find 


ea Catch a nat etek besE, vas we 
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t, Ww tk) 
& : s 
Bl atl inl eT s]-ty) bus d 
-5 (" de (ager? Fei wmitmes Lae] PP # By COU” |seas 
3 
& 
a f atf 6,84 dV = 0. ((2.16) 
° y’ 


where n is the unt outward normal to S aed Cy; ] 4s the jemp 1 € 6 aer0s5 5. 


Frou (12.16), we have the Euler equations ww V: 


Teh e CUM Ta gasgd + 6 Cag gs a 


es ee - 4 + Ee? =0 (12.17) 


and, 10 vi‘: 


Pit = 03 (12.18) 


(12.19) 
n(-60 ei] + PP+ RP + gu’) = 0. 


As may be seen from (12.16), boundury Conditions alternate to cach of 


(ang) are the syecitichtion of ui P© aud p, respectively. 


tit Constitutive Equations. 
We consider the constitutive equatcons for a cry shel with Macl-type 


lattice structure, whieh Falls under cube powit group mam, the qe evators for 
whith are given th (70). Owing te rhe turst gencrator in (hilo) ( He centre~ 
symmetry generator) there can be no coef ficreuts of odd rank and no prosucts 
of symme free aud ante symmetne variables in the qe functoon W. This 
reduces the interactions fo these between variables connected by full lines i0 Fig. 
1@) wheveas the elashed line segments deuote excluded wateructions . The 


Corres ponding energy functrow ts 


Se TEES BAP TEE EE Mer OUTS FECTS EL WATER EE OREO has RAD 


S ors oa SERS vy » , x, 
ats 3 oe 7 ae ms bs at EE SNe Ed 
a ie Pe eee ea Le: aes es ae. Oe ee ye eee 


PEELS EDLC, Aerts ene in 
{ 


Re gel ai tote A ic i 
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= 1 t pil pw p, a) ae e)_ A) () rah 
e Wee az 44 (a; ‘ae PF . bine P ic Pa * Ccjee Si Spa + 2d ine Pou ka ) 
#9) © te * 
+ Z (as wut P t2d** wi tii) * bate +c cog; We 
0 —&®) 

: +Z(b (0 Se F (2,20) 
where Pee) aud es clenote the symmetric aud autisymmetrce parts of Bese , 
spectively. We note Het, ta general, 

KA. Ale kA Ak kA Ak 
Air? ee y Deine > berks » Che = chee » Aejee ® Ties (2.29 
and, from (7.10), 

KA. KA 
ae: it a 

kA KA 
bing = 0 6 ka F be oe, Bie + bg (58 abet Sa Syn) br Gig be Sbinly e 

KA. KA WA : 
Cisee © Sa * Sus t Cae ae it Uk)» 

. A KA 
dies =A Sie + Un Oe bec + dan (Sig dig * Sey Sip), 


where oF (er Jeske) is unity when its indices are alike and zero Other wises 


bX = ‘A 


KA KA KA KA KA kA WA kA {RA ka 
a ~ 4, 2b, Shy Hh 2g, Ld edy = ey - 2g, 3 (12,23) 


a ua P . ; ee he 
for eee Chie and dan) the abbreviated notution for pars of widues, (jy or 


ke, has been used — asin Avtrele Ul. C+ mee 4 bec observed tat material 


constunts with superscripts I] or 22 denote interactron: within one of the 


two Component contimuu whereas Constants with superscripts 12 or 21 dente 


teractions between tre twa component Continua - 
From (12.12) and (12.20) - (2,22), we fiad phe Coustitutewe eg uations 
Ts MG ECG Le See Cit Ge Spt Rigg Se) (ds va en ta Gj on Techy Pe) 
. : Ha re aus, 
es sb G ZU ag Pe the diiPage tZbyg Rit, + Zbry Ay) 
nate § ECA Sys) oi) ea 5; 5 adi 52), 


i? Za* oi tau, 


2 pa" an r2c** 


(é.24) 


=~ 
: 


we . 


a tg ah ie oe <7 
LENE HG RE a ip te ge OL Gag a oe er ge REESE Mm Ae Appr, ag a 
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We shall assume thel, in tHe inital state, te. when 


ta} te) (a) + 
Pea Pecans Cigar @) se 
C4 Ne Me MG OF, Eb, 


a 
B 
ee 


theve t's “oOo 


res lfent . . 
. at force ocross aay surface — exferiwr or watereor . Thafes 


«) ¢ 3) 
t; rt; = O, 


Then, from the first of (12,19) 


: 7 7) 
3 Hele, Ay ty =0; 


SE ARTEL RPE 
Seca enttnctie ras Heit he tain ate eee tane co Gas es hk eu 


sepscaen 


SHS Tg 


- 206 Pa, 
core’ =O, (225) 


Accordingly, from ((2-22) end (1225) th the instal shete : 


Succes 
See na nS ak aS err 
. 


4 a | -W > 
% Celis Ty = 9, tat Ti? = 0, 
2 but there is a self -equili brated 
i EW. EH 
E; 4 6 be 
a ' : 
ES i i each component continuum. 
RS H 
i 
j 
: 


4 : iV. Surface Energy of Deformation anc Polarization, 


x 7 
4 Ln a state of eguilibraim, the total energy in Ve i9y- Peoitg C013) 

| a 

q We: be 4 
a fw dv + 2 i 4.9 dV, (12.26) 
3 Weean write W in the ian 
rs { 2 &. a Col (Wo) (e) aw é ¢ 

E We ET Se AME RA) TT IUt eT oad eZ ORY nos 5) 
% | Then, by the some procedure os Mat employed 14 arriving at (12.16), 

E ra >. wo , Ke ® 7 * os ; 

j | Gwin sz OTe, ae (ePeeap 

+ a Ca) |e 7 #7, ot) 
FE | % Semalrge Cat ret]u) r Enemas +3 E fini (oR actus Jas, (2.27) 

: f 
Oo 
t ; 
— ee 

-.. Be co a5 GE “Re ae reece eye tere LT Or ne RRR saree ucaiaeaies Fe Fees capt tes ae ve 


Ba Ee ie ie ete ener 8 Fe A baa a Tah ile 


BN A SEP LIN EE OE RE EI ANSE TOSS EY ELT UTE OE 


a Nac Et De ae ae ks SR RAE EEN EGA cake calle eo ae be aes chon eae a calccy | 
Seat 


ie 


17 ; 

| 
Also, 
eas : j 
Le tegidve [oa Gish V# [ nelaile ds. (12.28) 
3 
Laserting (12.27) and (12.20) tn (12.26) aud applying the epan tons of eguilibraing 
i t.@- (12.47) wth che > 0, in the absence of external forces ‘cs GC, +0) we fiad j 
! : ; 
: { 
i : te! ( : : 
2m =-El Cog, ay re Rt (Portas de )Q]aV | 
q 
i ) : 
; t Ef mili +(-1)* Tet Ju; ren P™ rel ¢; lyjas r Ef melo" e, WMS. 
i 
a” § 
: Now 
wy 5 i 
fabaedye-( chav ef me Peds, ) 
, 
fy (Ure, a CAV=f, Megeate AS. : 

Heace 


2k ip n{ LTP COT UM rE e+ (elec ]- PM -@ ute}ds 
+E f. nz (bt P re*u* dS, (12.29) 


From (12.19) the first inte gral in (12.29) vanishes ¢f the boundary ts entirely free . 
= = ¥ 
awe Z fn, (eRe rciu* ds . 


According! , the Surface energy of deformation aud olanzateon ruurk evea ts 
“4°71 ie P 2 Ve : 


Ss. 
Wasa le PY eb Pm ecu" s. (2.30) 
Thi's resu(t ts to b- Lom pavect w tt (10.9). : 
v. Wave Motion. 


The eguations of motton iH terms of ge Pr"! aus @ are ob farmed 
by substituting the constitutive equations (12,24) inte the Field equa tous (12.17) 
wtte the result for uch tk, 
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ZS) KA A) re) 
pa at Beica Ce bi + Cy U Bis * ge (Up ti rail 


+ * [a Sits Pe +k p& deg (PY ithe =) 


oi * oy 

eZ Carga re aft ui | 
(] e oo(x) 

* (-1\* Lo" p? rattp” rat* (yu s " +9.) -9, 65 res 2 pie 


‘ (12.31) 
Fd ine os due us tg + dg ( 4; RU YEON (~t) d* (ur “a my 


+*Ze"S, aie Pye tein Poss ros’ (? a} sae (94 - p. 


oxy) (Pai ity jit UY 


te W_ @ a QQ) on es ‘ 
a, FP ay, B® - (um ui’) w+ FP =O, 


no 


10] {23 (2)_ ayo ; 
~ $F + Ri * hi + Gul 4; ui) = 0, ‘ 


The essential features of wave motwn are revealed by au Cramisation of 


plane waves in, Say, the x, , or (190), direction + 


te) aw e Ua -wtl ps Bie eee) ee Cel (te-wrt) (12.32) . 


3 
and (, Eps 0. With these functious, the five epuartions (12.31) reduce to thae 
for longitudinal waves if j=l and transverse waves if j22 or 3. Ths: 


Longitedinal (j21) 


17) Hp) #258) 887 (2) Ut) (ee) 
Cu in tn din thy Py tan Py take eo be ea™ (uy, -e() +9, By ga ie 
az (a () Hot) 4254) (88/7 (y ,,G) rey cal 
Cu Mun Cy Gy edn Po tds P ra -g Re -a P, -a (u: <4, }- Ke Er = Ff 77 > 
cu) 2) Wp) 3 pO ole td, rr 
di Mya rd Ayu +b, fu ’ bie Ru o aii, 9 -a\, F a (uy; U, } €, = 0, (12,33) 


at edit Uy 46H Pe) eb Pa — hn Boma PY AM (U-a0") -@, =0, 
aN + ane Bel, ~ yi" ) - Chu = 0. 


J 
_Transuerse (2 2) 
) 
(cl, a ee (ag ecu ra CM -d*') fei + (as ~d??) pe Ae 


tpt : a” er.) ig | Us ay vee Ww te 


3 a 
2 q c at) (a 12 at) pW) a2 
(ye ch) agi + Cog ~ Duy + Grd Bn . Ge, +d") Bay ay 
@, ¢¥/,, (2 WW) . , (ab esta) 
~q*! pi -a*? A - (u; -y.") = pt ‘ 


(12.34) 
(dd) vgn + i rd yy + (byy eb) B Ry + (ogy » bn) Pan 


tp (t) at Gy) . 
- ay ‘ee =i ie Cus Yas) = 0, 


ou ays sas 


¢ 
4 


5 
7 


oN 
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a jez) az 432) ai ar 2\ pu 42 @ 
ido" JUgy + yy 4 )Uay + Coy, tbr) Pay + (oar bi) Pen 


- ai! po. qn 2 - a (ur roa )= 0. 


I+ may be noted that the Maxwell self-Field ts coupled to he mechenica ( 
clis placement in the longitudial waves, but not in the transverse waves . As 
ts Shownin the Following Article, if the magne tee Freld 15 tekan wte 
gecount, te! longitudinal wave ts unchanged but tere 1s ful! coupling 
ot +he Ela ste and elect romagnetec fields in the transverse wave - 


Upon subs tituting (12,32) im (12.33) aud elimimating sige B,“ and G, 


Bags 


we find the dispersion velatron for lougrtudiaal waves 


ey 
a aN fa 


* 


Le eNO ACORN MEARE Fe ate vhccanpied i aktale abled 


AL= 0, (12.29: 


~ 
ov 


pies are AL eae 


. 
Laeyenie res 


where AO, is the determiant with elements 


| aie pew a**- gre cuss Aa, -o"'~8, 6," dig 2452, 
| 4n- arg es 7 aig = Bai, Bry?” -9, é; -divettay., 
Agea" + E66e ~ die" =O, Ay? -a, - & ~%, &*, (2, 36) 
Ai, = a +96 -div $* FO, A34 = - di - 6! -B 8* = Bas, 
| Bag? pow -a- ai Sy Dee? Hi -G HE 
i This is @ quartic eg eaten in &*. Of the fosr branches, two are rtul: one 
. acoustical Gud wwe optical a> mtu be verified by showing that one real 
| oranch passes Harough $=0, wi0 wud one +h rough $26, wel. We tnd, 
i” fact, 
eee te (eit ect ezeit}e® - (pp) wD, , 2,37) : 
; wheve ; 
| aivel gisele aye. ; 
Do= | ane ai re; a’*+ 9,6," |. (12.38) 3 
| a's a6 a+ ts! a*+gs é,' ; 
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. ae FRAY PORE a WERE EBSD SEN Eee ee a Ue es MQ gio aD SGA Sia stir A Be Ai aye sb Betty cee 
Rime ehye Mtr ea He 


a einen ee ee 
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Hence, at long wave lengts, the free uEREN of He longitudinal acousteé wave tS 
wine EL (en eC e2e, peep], (12.39) 


It may be seen from (12,39) thet the usual low Frequency extensione/ stift- 


Ves>, Cu, ¢5 given by 


A 
A 
| 


Cae Ge roe e 2c)! (12.40) 


; 2 
ce. the sum of He exteasional st, Ffaesses | cand cy of the twy Component 
Continua augmented by the satercomponent extewsiona! stir fuess 2c 


The long wave lait ef the longi tucine ( Optical branch 13 obfemed 


from 


lim A. zw fp et [ (ai ré. )( a2? +6) -(ae +6) J (p"+2")Q3, (1z.#) 


Here, the limiting Frequency of the longi tucine optical brach 1s 


ible toby Wig nt ade te eS ELE a DAR mI Pe BOIS TU CARL ORUE WE 


wr = {ea /BlUeare lean e)- (ir are} (i242) 


where 


p= e/ (pe aap (12.43) 


We. b is he reduced density, 
Ina similar we WV; we find, from (12.34), the Aspersion relation 


for trausverse Waves in the (100) direction: 


° 

P 3 

3 4,=0, (2.44) 
a wheve Ly ts the determinant with elements 

a ) 2 i s 

7 ae cs 2a) "= (C49 a Ays= “a ' (dy ed*')e* 452, 

a Ane O° ~ (Que cE = Aa, Aig =~" - (dag rd”*)8? =a,, | 

F Bye at - (dy -d°)87= Ay, O33 = a, ~ (oe, + Bh) 8", (12.45-) 
‘ 44° Cale das ~d**)§* * Bai, 43, > - ay * (ye +b. 7872 D4 

4 Daz * Pat ahh (Cae )5, Bag ® Oy -(ey rb) 8", 


Wer, Twas 283} 7-7 


, SWAY LAS LES Ltt SY PTTL AD ATRL E EIS RTI TIRE OI Tt RR BL FE 
<r PTET BC TALE SLORY FEN FER ee GSP EMA NG SE PSB AA SS OLEATE TEIN Re HE AEE we : 5 


&l 


Again, thes (8 a quartic in go, with two real branches. For the long weve &- 


Lae Ok, Oe 


havier of the transverse atoustic branch : 


* b bet 2: 2 q) Cs) 12 
im, 4,= Geek (ened t 2eg2 )s = (p rp Yeas JD, ’ cd 


Fe rd malted ental hl Ste ech tt Ni RRA NS ed 


where 
ay, a; Q"' 
D=la, ae al, (12.47) - 
at’ ot a” : 


Itence, there is a transverse acoustee branch wrth Minetvg, low Freguency 


keha vier 


se PR a 


=8 N (ayee r2 coe) (peo) : ; (12.48) 
The long ware limit of the treusve -s€ optica/ branch «3 found From 
fen, Aye opp w'( Qian dn )~ (pp )D] (12.49) — 


frove whieh he lisa shed Fregueucy of the fausverse optical brauch ts 
~ ‘s O 1/2 
Wry = [0/5 (aa; ay +a, "any (42 70) 


vi, Cne~ dimensional NaCl-4 pe Lattce of Shell-model Atoms, 
The purpose of thes secheon 15 to show thal the POPE POTEET 


equations of longi tadenal moftom (12,33) are the louy ware lint of the 
ditterende equatcans of matron of a cue-dimecusconal We Cl- tye lat tree 


of shelf mode! atoms. Such a lettne 1s couvenrently reyre seated by two 


ties SS 


ienes of alternating atoms (postive ancl N&getive fous, 19 the case ofaa' 


raph, ote 


: a (ome cristal) with one atom of each type aetcach leftice sche, as Shown 


I Freq. 9 where the two tuyes of atom are rdentifred by the arg, ts f 


Skt iy bit tenia bed eal 


oR as 
- 


cud 2, Meerest neighbor interactions between unlike (eed jaceat | atoms in Cach 


hae and next neertat necghbor jnteracttons between like atoms co the two ieee 
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Same site are cousidered» With each of p aud g teking on’ the identificatrons | 


or Z, the force constants of Ke interactions are denoted by dp for He wtrar 
atomic core —~ shell lirferactoous aud bevmecty for the care — core, core — shell 
aud shcll— shall piers tnferacaens vespecrively, ce ibe afous for p= 
4 aud unkke atoms for p#Q. We note that Boy > Pepe aud by = bp bat fig # ir: 
The two lines of atow's are takeu te be perelle( to the x, axes with tHe atom 
sifes at x na, where nS @ positive or neguhes ui teger. The Stsylace— 
ments of the cores Gud shells of the avon. wf Kona ave deubled by 
ut aud se, respectively, where eel, fo destqna:e te afous of types {- 
aud 2. Cr for the positive t ton aut 2 for the ugative fon, at the case ef an 
‘onde crystal). 
The equation of wotion of the n™ atom of type | 15 ob fumed 
by equating rts inertia force to the suum of the forces on its core aud shel/ 
as by He coves aud shells of the tyro nearest necg kbar (unlike) atous, 
He two next nearest necgkbor Chike) atoms and ihe Maxwell, electric setf-fiedl 
at K,2naS . 
Bur (unr “Un ) Yin (Spor~ Sut Une “Mn d+ dy (sy2, ~ 50°) 
~ Bin (on eas) = palettes # Sy ) — Syl ses, ) 
+ parla ce") ee (sm a") + so -56)) + Galspet ~5e ) 


~ pilua- ue) Teal dys 5pe) > pals er ~Gia(sn Se) “ee, emit . (asl) 


Similarly, fov the shell alone, of the type | atom : 


aa) latest) rd late) lau) By) 
+ prs(Unes 3a” Gal ser= 85°) “Yar (Be cge,) = Sig (SSH 4 BE =O, (12082) 


where % ts the decfronic charge and the mass of the shell ts neglected , 
We adopt the following detintius of electronc polertzation deanties aud 
mass dense teds + 


| p, (50° un Mela’, po =m, la? > keL2 | | (12.53) 


a ee ee tee meme em ert re nah neem es nites Ree ey ae eee eee 
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and also make the Following (dentifecateous : 


a, = (a, tery tl yar §.)a°9", ie gia a'/e = fa? 3 ae 
by = oi oa, a = aa? = (Bia Ya ty tS, la”, 
4 jecz. (12,54) 
Cu = (6,, tam + b,)a ’ =a af= (rar dala’g™, 
di = (yi + &,)a'9", : =~ 70°72 = (ry +6 )a'g" 


Then (12.51) @ud (12.52) can be witten as 

Se WE AT A PY dare PY. #? pO sg ra Cus us')-4,E, = pM,” : (12.55) 
di au” rho Aue bi i dues b? cee 2 qa" p -@ i pel. au. uy”) rE : 0, (1256) 
where A” is, agen, the second centml difference divided by Pa 


"hh the same way, for the tyye 2 atom, 


OU! 6 Ue aR raat PY ghpo_gitp®. gM (y8_yil) te,E,=p%ay , (57) 


dy Buy rd 0 yt by, BPA + APO ai P.O -g Pg uu) + E, = 0, (1258) 


where, in addrton te ((2.9°4), we identify 


ay = (a, i 2 hat 2 Taté §2)a°¢*, b, = 24 *e", 


: (i259 
Cy = Brat 2p " b22)@ a di = (yn+Grle’g"' G: 


Also, as th Artecle 8, ‘Le Charge Puc tion aba be written as 
~E dd y, + [Pe PY re (unas) =o, (12.60) 


where Dd’ and D are the Taylor series expansions of the derwative O/Dk, 


in terms of Forward and backward differences, respectively, 
In te long wave limit 


ui u(x, ), pl» P(x), E.,-7 E,(x,).=-d¢/dx, 


(12.6 
AE > O*4(x,)/dx, , 0 > 2/dx, . are 


Then (12.66), (12.57), (12.56) and (12.60), inthat order, reduce fo the 
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five equations (12.33) 1% the ofder given. Since the continuum Oqua trons are the 


a long wave liunct of the lathce equations, the two dispersion velations arc ipse facto 
q the Sume in the long wave hat. 

It sholild be observed thet, strictly spealcing, the real segments of 

te optical branches tn the continuum ayprovimateen, have no valid role in 
Cases where the optical aud acoustic modes arc coupled ~ as, for examyle 
tw soluftons of the Gq tatrens for free boundary condi trons. The frequencies 
g ot the two Fu pes of mode must be tHe same, wher Hey are coupled | 
Dut, any Frequenc 4 ou a real segment of an optical branch 15 so high 
Hat the wave leugrs of the acoustee bruuch, at that fre wency, “s Too 
E Short ty be wcthin the range of validita of the continuum 4 yp proximation. 
a In such situatcons Hee 1s, i genera (, an upper lime # of frequency, 
E set 4 the wave length limitation of +he acoustec branches, tywhich 
rhe solution fs restricted, 
ie Howe ver, there ig an 1 por faut exception. At Frequencies cloe 
to Wyo or Wr, the wave length of the opt cal branch ts extremely long 
§ i” Lom parison with that of the acoustec branch, In that Case, the amp- 
: litude of the contributron of the Kcoustec pact tu the couyled mode 1s 
small in comparison mth the amyl tude contirbutton of the optecalpert ; so 
that the error tin tocluced by the thuccuracy of the acoustee branch 15 small 
4 In fact, the acoustic branch may be discarded , then, without seriously 


atfecting what remans — thus Steal Frying the equations of mutton. “he sct- 


uation ts analogous 46 Hat ia the theors of ttickness-shear aud Flexural 
vibrations of plates CoS). The theckncss shear branch corresponds to the 
optical bvanch and the Flexural bwanch corresponds to He acouste brunch, 
At Prquencee s near the cut-off Receeey ot the theckness- shear branch, 
He amplitude ratte of the Flexurul to the Hauckness- shear modes 1s smulf 
£65]. In vhat case, the Flevurel breach can be discarded entirely aud a 


Sim pler egue toon of motcon for tHe teckness-shear ode, alone, nay be 
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written [Les | 
In the long wave, low Fregue mey limit, we have i udditicn tg (2.61), 


2 WM 3 @ 
u, =u, , P, an : 


Then, with the notatcons 


~ 0. ,,00 ~ pw (2) a 
u,zu, =4 , 3 P, P, cA » (=P rot) 
a 22 12 “ 12 t2 
Ayr olan ray, r 2a, = b,, +4 (bi. +b, + 2 by, i; 
Cue Cut Green, d, 73 (d +d Sede ds) 


the sum of (12.65) and (12, s7), 


tre sum of (12.50) and U2.58) and, Fusatly, 


(12.609) heceme, respectively , 


CH dyy Fu Pra = Oa 


) 
ania +b, Fun > a,,P = ( 


u 
& 


> ae ai Pui 


These ave the equations (Ale) extibited ia Artete 7 for lenge tucevtue ( 


wuves in a simple, as opposed toa compouna, Aielectte contimuuu, with 


symmetry mom, wheu the polarizatcon gracient is taten wty account Thus 
the long wave hmit ef te equations ef the onc-dienstonal, Giatemce lattice 


yreld s the equations tev longi tudina waves in the Compound continuum 


when no restriction vs Placed en the trequene a7 whereas the long weve , 


low i eae lime t yields the €q,wutions se the simple continuum. 


tt very Vee Frequencies, the G@couste branch ts nearly exuct and, 


at zero Tits theteg (this @bsent. The qunlety et +e apprakwunation then cde - 


pends on the low or zero Fregutney values ot the ma qinary or complex 


branches a f He dis persion relation . These affect behavier 4h Or near 


surfaces. A typical example (5 given by the difference between ch ana cd, 
the magia avy wave lengths at zero Frequency ot the Imaginacg branches 


of the continuum and lattice dis perscon relations, in Article 8, ana the 


ae eo 


+d at ahs 
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role of this difference in the solutcons of the continuum aad lattices egua- ; 


tions For capacitance of Hin Films 1% Article 9, 


set rp eee teh aod BE CaP Teor s . i Korg iF PREM DX oa GPa fae ee 
Ba : ete: METALS REERU RAS RUS a we a Shake a= a TT Ae rts OTe INT AA aS a oy Ben se. isesa Sune 
FPP SEED & pies pe Oe AUG or eo oe TERNS : 

? 
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13. Coupled Elastic and Electromagnetic Fields iva 


Diatumee J Drelectric Continuum 


Equations governing coupled iceLanical and electromag netic fields an da 
tome , (once, optically isofrone erg stals were given by Huang [67,68] i terms 
ofa polars vation variable, the relutwe al ts placement of the tons aud the 
uguel variables of the electromagnetic Freld, These Cguatrens have been en - 
tended keke tate account the separate electron. polarizations wud. the 
separate displacements of phe tur tons, eke toned polarita tren, the two clectronce 
polar tation gradients aud the twy aispacement gradients , Asede From a 
mone detuled represey tation of He atomic and electronic wa teractions , there 
ave two man effects of the addrtonul considerations | first the equ Pans 
ave a pplecable to shorter wave lengths , ow ing fo the (ncluston of the dis- 
p cement und polarizatton gradients second, the acoustic brenches are 
: neluded ya the disperscon relations for p lune waves, re addetten to the 
optical brenches and the electromagnetee branch . 

As ta Article 12, each of He two ater penetrating continua repre set 
Ing the Statomec ery shel has (fs own As placement dey k=2, aud elec— 
troutc pelarre ateon dae The stored energy of deformafton auc polarizaton 
bs the same as before: u functron of He dividual Small shams , si, the 
tadiidual Polarizetions aud polarieatrou gradients, Pe aud ae and te rela- 


rive dis placemen ee as aud rotation, ww 


ey 
Q) 5 W) 5Q@) pW pa | *F 
whew * (58, Sia Pe Pr Be Pe He, HG) (13.0) 


4 


“Also, as he fore, the total polar;retcon, per unct aren P:, (s the tum of the 


electronc polarization s aud the loucce polarrec anes 
Pes ps" +p? + Geto : (13.2) 


CE the vane bles (4 (13.1), only P: and ue ave accounted fos iy Huang’s 


[e768] gua tions of motion. 
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When the Maxwell, electrc selt- field, Ey, ts quast- static iF satisttes 


7 


Six &nj = 2 (13,3) 


2 
where , as in (3.2), Otik 1s the unct alternating teusor. In the alec tromagnehe 
case, (13.3) ts replaced by 

Sik Eus + 8, =, (13.4) 
where, as in Article ll, the magn etec Flux density, 3; sakes tres 


Ay Sein bs > Ey of » 8 =O. (13.5) 


Ct 

e+ 
In addiizon to (13-4) and (03.5) the Field cautions (12.31) also hold, bit 
with A replaced by Ei 


Plane wares th the Ftoo] direction are expressed by 


uc" &) ec? @)  , Ge cy =o ty 5 iy 
> K. ; P. =Lie |, c; -Me@°, G.=N;e (13.6) 


Cs) Ce) 


where y= §x,-wt and Kz , Le, Me aud N; are constants. 


Tu the case of long: tudinal waves, 


is sus" = 0, ae ee =0, &,-F,+0, 43-7) 
u@akeel®, PMaLMe, Ermel cos) 


Theu, from (13.4) 8,205 Ge. there (5 no coupling with he magnetic 
Creld. Mowever the Maxwell electne seif-freld us coupled wet the 
displacement and polarization fields. The solutcon, therefore, (s the 
same as tat found in Article (2 for the case oX the guasc- state 
electrié Field, with dli’sperston relaton (12,35) as the result. 

Ot rhe two stinilar transverse waves, with ds placements 14 the 


divections of x, and x, , we chovse tHe former for examination + 


we 20, PM=P 20, E,>€,=8, ae (13.9) 
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Tun the case of the guasc- state electric field, the absence of (13.5) permitted Me 
E,, as wellas 6, , to be wero) but both must be non-eero when tHe full ae- 

: tromagnetic tguatcons are imposed. 

d To fkd a form of the disper seon determaent similar to (12.44), 1F 0s 

useful First to express M, and Ns tn tewms 0 © HM ant 1 Hrough tHe use 

f of (13.4) aud Cas): 

| Mae io Myge- eh [LP LD ee (KO WP), (3.1) 

where 

: = 6-2" a C13.12° 

i ee ew wwe ER TE oe ae 

4 . 

i then the results (13.1), along with (13.9) and (1310), are iaserted in Ke 

remaining eguatcons of motion, the (latter are sats tial i= 


E 


ay be 0, (13.13) 


where an is the determirant with clemeuts 


Au= Mw! a4 LC “s, day “0"- 4.6 deg LF 2A, 
Aa? ote gig" a (ce +c**)e2 = Ais bee? vo gee rd™)B A, 
A= a gee’ ~(dyy -d*)E? = 05, sy? -a,, ~6;'~ (by, HOS) EP, 

A= OG yey! (be S)$%2 O44 Oa? ar! ~ (egy bi) F2=dy, 


Ban= PM wi =a" ai es (ae - OES Bag = ~ Ou ~6;' ~ (eb + b3)82, 


The long wave, nigh Hreguenen behavior is obtamed from 


: hon as" > ofp oes" T (ai, + 6, ai* ite, Z ~(a)2 +6, y*] - ~ (pre) D3, (3, 14) 
; where D, fs again given by U2. 38). Thus, tle lihartens, Freguenc y (S : 
wr9 0 = 1D. /pllara'y( an sg &)- - (a;' re) Ty (1315) 


C.t., as Found by fHuang 161,683, the same as the luting Frequency of the 
longi tudihal optical breach . 


Be athe ele nadatinn band iva 6 
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. 
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The long wave, low Frequency behuwrar .s qiven by 


din, AS" = tea Ps Coat Cre td cya) e (pp) ax" ] Dy , C13.16) 
where 
Dy = (6 O w~ 8*DY/(€, 4, a - &*) (13.17) 


and D is agath given by (12.47), Thus , there are two branches, with 
limoting, low Frequency , long ware behaviors 


PEAT Re Ta Ry A ee aE RPE BaP ae ROB we gpl 


é t 
e Wen > § (D/ea. D.) | (13.4 8) 
wea = SL Gy Ge + 2c) (pre )] ae (343) 


| 
| 
} 
| 
j 
: 
! 
' 
I 
i 
| 
: | The first of these i's the long wave end of the elec tromagne tec branch, witt 
| slope equal to the joreduct of the velocity of Cleetromagnetic wares 
ina vacuum (6, as and the thder of yarncticn (bp, /0)"" ii the 
dielectne. The second ts the long wave end of the transverse acouste 
ranch : identrcal with that for the quast-statee electec Field (Artiele 12) 
and the purely elastee field (Article +). The long wave portioss of the 
| real branches of tHe dis persion relations for both the long) fudinal and treas- 
verse waves are (Mustreted wm Fig.20. The vertical scale ot the acoustic 
branches (6 expanded: The rl® of Herr slyes te tat of the elect maguefec 
brauck «3 of the erder of 10”. 
The loug wave bekavirs otf all the real hauches are imdeyeudent of 
| the coustarts bets, Cejka ane dejh0» ne. in depsndent of rhe polaritercon 
| aad dis yplacemen 7 gradients. Also, the loug wave mets of the a@couskee 
branches deyeud only ou whe mass deasities and rhe elastic str ¥fuerses. 
! Hence as fur as the line orig bekavrcaor at the long wave Cud s$ concerned, 
tte results, here, Con Yorm with Hnang’s, with the adde¥rou of Me acous- 
te branches. However, as may be seen From the determixanrs (12.35) 


aud C1363) | when the wave length Simrishes From in finite, the dis- 


Placement @ucd yo lari zattou gradients atfect the disperscon relatcon, @s 
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does the coupling with the elastic freld. Then the results dwerge trom 
Muang's 

Avery qeneral treatment, including crystals of arbitrary sym- 
metry, an arbitrary number of ions and electrons yer unit cell and 


non- Itneer, as well as [snear, in teractions, ts given by Lax and Welson 


[7e]. 


rue 0s 
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Attn: Mr. Garet Bornstein i 
Research & Development Div. a i 
Indian Hdad, Maryland 20640 


" Chief ‘of Naval Operations 


Department of the Navy 


. Washington, D.C. 20350 
Attn: Cone OP-07T 


Deep Submengende Systems 
Naval Ship Systems Command 
Code 34522! 

Department of the Navy 
Washington, D.C. 02360 


Attn: Chief Scientist 


Director, Aero Mechanics ; 
Naval Air Development Center: 
Johnsville i 

Warminster, Pennsylvania 18974. 


Naval Air Systems Command 
Department of the Navy 
Washington, D.C. 20360 
Attn: NAIR 320'Aero. & actaee: 
5320 Structures 
604 Technical:Library 
52031F Materials 
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Naval Facilities Engineering Command 


’ Department of the Navy 


Washington, D.C. 20360 

Attn: NFAC 03 Res. & Development 
04 Engineering & Design 

14114 * Technical Library 


Naval Ship Systems Command 

Dept. of the Navy : 

Washington, D.C. 20360 

‘Attn: NSHIP 03I'Ch. Scientists’ for R&D 
0342 Ship Mats’. & Structures 
2052 Technical Library ' 


Naval Ship Engineering Center , 
' Prince George, Plaza 
Hyattsville, Maryland 20782 
Attn: NSEC 6100 Ship Sys. Engr. & Des 
,6102C Computer-Aided Ship Des " 
vee 6105 Ship Protection 
6110 Ship Concept Design 
0120 Hull Div. } 
' 6120D Hull Div. 
6128 Surface Ship Struct. 
6129 Submarine Struct. 


Naval Ordnance Systems Command 


Dept. of the Navy : 


Washington, D.C. 20360 
Attn: NORD 03 Res & Technology 
035 Weapons Dynamics 
9132 Tech. Library ‘ 


Engineering Department 
U.,S. Naval Academy 
Annapolis, Maryland 21402 
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Commander WADD 
Wright - -Patterson Air Force Base 
Dayton,, ‘Ohio 45433 ' 
Attn: Code WWRMDD AFFDL (EDDS) 
Structures Division 
A®LC (MCEEA) 
‘Code WWRC AFML (MAAM) 


Commander ; 
Chief, Applied Mechanics Group 
U.S. Air Force Inst. of Techn. 
Wright -Patterson Air Force Base 
‘Dayton, Ohio 45433 


. ‘ 1 
Chief, Civil Engineering Braach 
WLRC, Research Division 
Air Force Weapons :Laboratory 
Kirtland AFB, New Mexico 87117 


Air Force Office of Scientific Res. . 
1400 Wilson Blvd. \ 

Arlington, Virginia 22209 

Attn: Mechs. Div. 


NASA 

Structures Rescarch:Division 
Naticnal Acronautics & Space Admin. 
Langley Research Center , 

Langhley Station : 

Hampton, Virginia 23365 

Attn: Mr. R. R. Heldenfels, Chief 


. 1 ! 
National Aeronautic & Spate Admin. 
Associate Administrator for Advanced 
Research & Technology 
Washington, D.C. 20546 


Scientific & Tech. Info Facility 
NASA Representative (S-AK/DL) 
P, O, Box 5700 — 

Bethesda, Maryland 20014 
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Other Government Activities 


Technical Director 


Marine Corps Develop.nent & Educ, Command 


Quantico, Virginia 22134 


Director 

National Bureau of Standards 
Washington, D.C. 20234 

Attn: Mr. B. L. Wilson EM 219 


National Science Foundation 
Engineering Division 
Washington, D.C. 20550 


Director 

STBS 

Defense Atomic Support Agency 
Washington, D.C. 20350 


Commander Field Command 
Defense Atomic Support Agency 
Sandia Base 

Albuquerque, New Mexico 87115 


Chief, Defense Atomic Support Agency 
blast & Shock Division 

The Pentagon 

Washington, N.C. 20301 


Director Defense Rescarch & Emer. 
Technical Library 

Room 3C -128 

The Pentagon 

Washington, D.C. 20301 


Chief Airframe & Equipment 
FS-120 

Office of Flight Standards 
Federal Aviation Agency 
Washington, D.C. 20553 


Chief of Research and Development 
Maritime Administration 
Washington, D.C. 20235 


Mr. Milton Shaw, Director 

Div. of Reactor Develo. & Technology 
Atomic Energy Commission 
Germantown, Maryland 20767 


Ship Hull Research Committee 
National Research Council 
National Academy of Sciences 
2101 Constitution Avenue 
Washington, D.C. 20418 

Attn: Mr. A.R. Lytle 


PART 2- CONTRACTORS AND OTHER 
TECHNICAL COLLABORATORS 


Universities 


Professor J. R. Rice 

Division of Engineering 

Brown University 

Providence, Rhode Island 02912 


Mr. J. Tinsley Oden 
Dept. of Engr. Mechs. 
University of Alabama 
Huntsville, Alabama 35804 


Professor R. S. Rivlin 

Center for the Application of 
Mathematics 

Lehigh University 

Bethlehem, Pennsylvania 18015 


Professor Julius Miklowitz 

Division of Engr. & Applied Sciences 
California Institute of Technology 
Pasadena, California 91109 


Professor George Sih 
Department of Mechanics 
Lehigh University 

Bethlehem, Pennsylvania 18015 
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Dr. Harold Liebowitz, Dean 
School of Engrg. & Applied Science 
George Washington University 

725 23rd Street 

Washington, D.C. 20006 


Professor Eli Sternberg 

Division of Engr. & Applied Sciences 
California Institute of Technology 
Pasadena, California 91109 


Professor Burt Paul 

University of Pennsylvania 

Towne Schoolof Civil and Mechanic: 1 
Engineering 

Room 113 Towne Building 

220S. 33rd Street 

Philadelphia, Pennsylvania 19104 


Professor S. B. Dong 
University of California 
Department of Mechanics 

Los Angeles, California 9C024 


Professor Paul M. Naghdi 
Div. of Applied Mechanics 
Etcheverry Hall 
University of California 
Berkeley, California 94720 


Professor W. Nachbar 
University of California 


Department of Aerospace & Mech. Engrg. 


La Jolla, California 92037 


Professor J. Baltrukonis 

Mechanics Division 

The Catholic University of America 
Washington, D.C. 20017 


Professor A. J. Durelli 
Mechanics Division 

The Catholic Univ, of America 
Washington, D.C. 20017 
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Professor H. I. Bleich 
Department of Civil Engineering 
Columbia University 
Amsterdam Ave. & !20th Strect 
New York, New York 10027 


Professor R. D. Mindlin 
Department of Civil Enginecring 
Columbia University 

S. W. Mudd Bidg. 

New York, New York 10027 


Professor A. M. Freudenthal 
George Washington University 


School of Engineering & Appiied Science 


Washington, D.C. 20006 


Professor B. A. Boley 

Department of Thco. & Appl. Mech. 
Cornell University 

Ithaca, New York 14850 


Professor P. G. Hodge 
Department of Mechanics 
Illinois Institute of Technology 
Chicago, Illinois 60616 


Dr. D. C. Drucker 
Dean of Engineering 
University of Illinois 
Urbana, Illinois 61801 


Professor N. M. Newmark 
Department of Civil Engineering 
University of linois 

Urbana, Illinois 61801 


Professor James Mar 

Massachusetts Inst. of Technology 
Room 33-318 

Dept. of Aerospace & Astro. 

77 Massachusetts Avenue 

Cambridge, Massachusetts 02139 


Library (Code 0384) 
U.S. Naval Postgraduate School 
Monterey, California 93940 


Sadie nena cotta at ae 


SAH 


red 2 


RUSE TT TE Fig PEER IPP SIN SERA Se Pas Bn? FERRIER ESE UATE US PEPER 
ee 


RK Rm ae ie meme 


Dr. Francis Cozzarelli 

Div. of Interdisciplinary 
Studies and Research 

School of Engineering 

State University of New York 

Buffalo, New York 14214 


Professor R. A. Douglas 

Dept. of Engineering Mechanics 
North Carolina State University 
Raleigh, North Carolina 27607 


Dr. George Herrmann 

Stanford University 

Department of Applied Mechanics 
Stanford, California 94305 


Professor J. D. Achenbach 
Technological Institute 
Northwestern University 
Evanston, Illinois 60201 


Director, Ordnance Research Lab. 
Pennsylvania State University 

P. O. Box 36 

State College, Pennsylvania 16801 


Professor ]. Kempner 


Dept. of Aero. Engrg. & Appl. Mech. 


Polytechnic Institute of Brooklyn 
333 Jay Street 
Brooklyn, New York 11201 


Professor J. Klosner 
Polytechnic Institute of Brooklyn 
333 Jay Street 

Brooklyn, New York 11201 


Professor A. C. Eringen 

Dept. of Aerospace & Mech. Sciences 
Princeton University 

Princeton, New Jersey 08540 


Prof. S. L. Koh 

School of Aero, Astro & Eng. Sci. 
Purdue University 

Lafayette, Indiana 47907 
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Professor R. A. Schapery 


Civil Engineering Department 


Texas A & M University 


College Station, Texas 77840 


Professor E. H. Lee 

Div. of Engrg. Mechanics 
Stanford University 
Stanford, California 94305 


Dr. Nicholas J. Hoff 
Dept. of Aero. & Astro. 
Stanford University 
Stanford, California 94305 


Professor Max Anliker 
Dept. of Aero. & Astro. 
Stanford University 
Stanford, California 94305 


Professor Chi-Chang Chao 


Div. of Engincering Mechanics 


Stanford University 
Stanford, California 94305 


Professor H. W. Liu 


Department of Chemical Engr. & Metal. 


Syracuse University 
Syracuse, New York 13210 


Professor S. Bodner 


Technion R. & D. Foundation 


Haifa, Israel 


Dr. S. Dhawan, Director 
Indian Institute of Science 
Bangalore, India 


Professor Tsuvoshi Hayashi 
Department of Aeronautics 
Faculty of Engineering 
University of Tokyo 
BUNKYO-KU 

Tokyo, Japan 
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Professor J. E. Fitzgerald, Ch. 
Department of Civil Engineering 
University of Utah 

Salt Lake City, Utah 84112 


Professor R. J. H. Bollard 
Chairman, Aeronautical Engr. Dept. 
207 Guggenheim Hall 

University of Washington 

Seattle, Washington 98105 


Professor Albert S. Kobayashi 
Department of Mechanical Engineering 
University of Washington 

Seattle, Washington 98105 


Professor G. R. Irwin 

Dept. of Mechanical Engrg. 
Lehigh University 

Bethlehem, Pennsylvania 18015 


Dr. Daniel Frederick 
Department of Engrg. Mechs. 
Virginia Polytechnic Inst. 
Blacksburg, Virginia 24061 


Professor Lambert Tall 
Lehigh University 
Department of Civil Engrg. 
Bethlehem, Pennsylvania 18015 


Professor M. P. Wnuk 

So. Dakota State University 
Department of Mechanical Engrg. 
Brookings, South Dakota 57006 


Professor Norman Jones 
Massachusetts Institute of Technology 


Dept. of Naval Architecture & Marine Engrg. 


Cambridge, Massachusettz 02139 


Professor Pedro V.Marcal 
Brown University 

Division of Engineering 
Providence, Rhode Island 02912 


Professor Wexner Goldsmith 
Department of Mechanical Engineering 
Div. of Applied Mechanics 

University of California 

Berkeley, California 94720 
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Professor R. B. Testa 
Dept. of Civil Engineering 
Columbia University 

S. W. Mudd Bldg. 

New York, New York 10027 


Dr. Y. Weitsman 

Dept. of Engrg. Sciences 
‘fel-Aviv University 
Ramat-Aviv 

Tel-Aviv, Israel 


Professor W. J. Pilkey 

Dept. of Aerospace Engrg. 
University of Virginia 
Charlottesville, Virginia 22903 


Professor W. Prager 

Division of Engineering 

Brown University 

Providence, Rhode Island02912 


Industry and Research Institutes 


Mr. Carl E. Hartbower 
Dept. 4620, Bldg. 2019 A2 
Aerojet -General Corp. 

P.O. Box 1947 

Sacramento, California 95809 


Library Services Department 
Report Section, Bidg. 14-14 
Argonne National Laboratory 
9400 S. Cass Avenue 
Argonne, Illinois 60440 


Dr. F. R. Schwarzl 
Central Laboratory T.N.O, 
Schoennakerstraat 97 
Delft, The Netherlands 


Dz. Wendt 
Valley Forge Space Technology Cen. 
General Electric Co. 

Valley Forge, Pennsylvania 10481 
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Library Newport News Shipbuilding Dr. Kevinj. Forsberg, Head = 4 
& Dry Dock Company Solid Mechanics 
¢ Newport News, Virginia 23607 Orgn 52-20 , Bldg. 205 
a Lockheed Palo Alto Research Lab. 
. Director Palo Alto, California 94302 : 
Ship Research Institute : 
Ministry of Transportation Dr. E. M.Q. Roren ; 
700, SHINKAWA Head, Research Department ; 
Mitaka, Tokyo, Japan Jet Norske Veritas : 
Post Box 6060 F 
Dr. H. N. Abramson Oslo, Norway ; 
Southwest Research Institute : 
8500 Culebra Road Dr. Andrew F. Conn 
San Antonio, Texas 78206 Hudronautics, Incorporated B 
Pindell School Road 3 
Dr. R. C. DeHart Howard County ; 
Southwest Research Institute Laurel, Maryland 20810 i 
8500 Culebra Road ‘ 
San Antonio, Texas 78206 j 
Mr. Roger Weiss 
call High Temp. Structures & Materials 
Applied Physics Lab. 
ad 8621 Georgia Avenue 


Silver Spring, Maryland 20910 


Mr. E. C. Francis, Head 
Mcch. Props. Eval. 

United Technology Center 
Suimyvale, California 94088 
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Mr. C.N. Robinson 

Atlantic Research Corp. 
Shirley Highway at Edsall Road 
Alexandria, Virginia 22314 


Mr. P.C. Durup 
Aeromechanics Dept., 74-43 
Lockheed -California 
Barbank, California 91503 


Mr. D. Wilson 
i Litton Systems, Inc. 

| AMTD, Dept. 400 

El Segundo 


9920 West Jefferson Blvd. 
Culver City, California 90230 
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